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ABSTRACT 


A  numerical  study  was  conducted  to  predict  the 
combined  convective  and  radiative  heat  transfer  rates  on 
the  walls  of  a  small  aspect  ratio  cylinder  representative 
of  the  scaled  model  of  a  rocket  engine  combustion  chamber. 
A  high-temperature,  high-pressure  environment  was  simulated 
in  the  cylinder,  with  gas  velocities  at  low  subsonic  levels 
typical  of  the  conditions  leading  to  the  entrance  of  the 
nozzle  section  of  a  rocket  engine.  The  composition  of  the 
gases  in  the  cylinder  was  determined  from  the  TEP  program 
for  the  burning  of  rocket  fuel  at  typical  values  of  the  0/F 
ratio.  The  thrust  of  the  study  was  to  determine  the 
radiative  contribution  to  the  heat  transfer  rate  from  the 
hot  participating  chamber  gases  to  the  cooler  wall.  The 
calculations  were  carried  out  using  the  commercial  CFD 
package  CFDACE,  and  were  first  benchmarked  against  known 
results  in  the  literature  for  the  simpler  case  of  gray 
chamber  walls  and  a  gray  participating  medium.  The  non-gray 
computations  were  subsequently  carried  out  using  gas 
absorption  coefficient  values  obtained  from  the  exponential 
wide  band  model  with  the  help  of  the  fire-modeling  program, 
RADCAL.  The  effect  of  different  chamber  wall  temperatures 
and  gas  compositions  were  examined.  The  main  findings  of 
the  study  are  that  the  radiative  contributions  at  the  high 
gas  temperatures  being  considered  are  comparable  to  the 
convective  values,  and  strongly  spectral  in  nature.  In 
addition  these  radiative  fluxes  were  found  to  be  least 
sensitive  to  the  wall  temperature  and  chamber  pressure  in 
the  range  considered.  Furthermore,  this  radiative 
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contribution  reaches  a  maximum  at  a  unique  optimal  optical 
thickness  of  the  gas  that  lies  within  the  extremes  of  the 
optically  thin  and  thick  limiting  cases. 
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I. 


INTRODUCTION 


The  heat  transfer  behavior  in  many  high  temperature 
applications  involves  a  combination  of  radiative  and 
convective  heat  transfer  modes.  Analyses  of  such 
applications  are  often  complicated  by  irregular  geometry, 
turbulence,  and  combustion. 

Radiation  within  a  medium  containing  products  of 
combustion  is  dependent  upon  the  temperature  and 
concentrations  through  the  entire  field.  The  energy  is 
distributed  across  the  infrared  spectrum  in  a  highly 
nonlinear  fashion.  This  greatly  complicates  modeling  of  the 
heat  transfer  within  such  high  temperature  environments. 

Gaseous  radiation  is  particularly  relevant  to  heat 
transfer  in  combustion  chambers  burning  carbon,  hydrogen, 
or  hydrocarbon  fuels.  Absorption  and  emission  by  carbon 
dioxide,  water  vapor  and  other  products  of  combustion  is  of 
particular  concern  due  to  their  strong  radiative 
participation . 

This  thesis  deals  with  the  characteristics  of  combined 
convective  and  radiative  heat  transfer  rates  on  the  walls 
of  a  small  aspect  ratio  cylinder  representative  of  the 
scaled  model  of  a  rocket  engine  combustion  chamber. 


1 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


2 


II.  BACKGROUND  &  MOTIVATION 


The  flow  and  heat  transfer  behavior  in  a  rocket  engine 
is  crucial  in  determining  its  performance  as  well  as  its 
thermal  signature. 


The  schematic  of  a  typical  rocket  engine  is  shown  in 
Figure  1 . 


Figure  1 . 


Rocket  Engine  Schematic 


There  are  three  main  sections  in  a  rocket  engine.  The 
first  is  the  combustion  chamber,  in  which  oxygen  as 
oxidizer  and  fuel,  i.e.  RP1,  are  ignited  and  a  high- 
temperature,  high  pressure  environment  is  created.  Gas 
velocities  are  at  low  subsonic  levels  typical  of  the 
conditions  leading  to  the  nozzle  section.  In  the  nozzle  the 
flow  is  accelerated  to  supersonic  velocities  and  is  ejected 
in  the  exhaust  plume. 


For 

this 

phase  of 

the 

research  program. 

this  thesis 

focuses 

only 

on  the 

heat 

transfer  behavior  in 

the 

combustion  chamber  and 

this 

thesis  provides 

modeling 

and 
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numerical  support  to  ongoing  Air  Force  Research  Laboratory 
sponsored  research  program. 
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III.  PROBLEM  DESCRIPTION 


A.  THE  MODEL  PROBLEM 

Detailed  schematic  of  scaled-down  model  of  cylindrical 
combustion  chamber  as  used  in  the  experimental  program  is 
shown  in  Figure  2 . 
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The  scaled  model  of  the  combustion  chamber  consists  of 


a  participating  (absorbing  and  emitting  but  non-scattering) 
gray/non-gray  gas  flowing  through  a  circular  duct  in  the 
entrance  developing  region  of  flow.  The  gas  is  assumed  to 
enter  the  duct  at  a  uniform  velocity  and  high  temperature 
and  pressure.  The  duct  wall  is  treated  as  isothermal  with 
temperatures  ranging  from  400  °K  to  800  °K  .  A  schematic 
diagram  of  the  model  is  shown  in  Figure  2.1.  The  inner 
diameter  of  the  duct  is  1.5*,  and  its  length  is  3*.  RP1  is 
the  fuel  being  burned,  while  oxygen  is  used  as  oxidizer. 
The  products  of  the  combustion  process  are  determined  based 
on  an  O/F  value  of  3  which  is  typical  of  the  chamber.  The 
actual  combustion  is  not  modeled  in  this  initial  study, 
although  its  net  effect  is  represented  as  a  well  mixed  core 
of  the  constituents  with  core  temperatures  ranging  from 
1500°K  to  3500°K  for  use  in  this  parametric  study.  The 
inlet  fluid  velocity  to  the  duct  is  taken  to  be  uniform  at 
a  value  corresponding  to  a  low  subsonic  Mach  number. 

B .  GAS  RADIATION 

As  per  the  laws  of  quantum  mechanics,  molecular  gases 
can  emit  or  absorb  photons  to  varying  degrees  at  an 
infinite  set  of  distinct  wavenumbers  or  frequencies.  No 
spectral  line  can  be  truly  monochromatic;  rather, 
absorption  or  emission  occurs  over  a  tiny  but  finite  range 
of  wavenumbers.  The  results  are  broadened  spectral  lines 
that  have  their  centers  (maxima)  at  the  wavenumber 
predicted  by  quantum  mechanics.  Three  most  important 
phenomena  that  cause  broadening  of  spectral  lines  are 


6 


collision  broadening,  natural  line  broadening,  and  Doppler 
broadening . 

A  single  spectral  line  at  a  spectral  position  is 
characterized  by  its  intensity  and  its  line  half-width. 
However,  a  vibration-rotation  band  has  many  closely  spaced 
spectral  lines  that  may  overlap  considerably.  The 
absorption  coefficient  for  the  entire  band  is  found  by 
adding  the  absorption  coefficients  for  single  lines. 

=  (2>1) 

j 

All  spectral  integrations  may  be  reduced  to  two  cases, 

1  Knhb)vdTl  and  £  I(b)r][\  -  exp(-£V/5)]^  (2.2) 


where  7(i)  denotes  that  either  Ib  or  I  can  occur. 


Since  the  local  radiation  intensity  is  due  to  emission 
from  all  locations  within  the  medium  and  the  bounding 
walls,  it  can  be  expected  to  be  smooth.  It  is  further 
smoothed  by  absorption  and  scattering.  As  a  result,  the 
previous  equation  simplifies  to 


1  M7+A17/2 

[At] 


K/T] 


drj 


(2.3) 


and 


(2.4) 


In  order  to  find  narrow  band  values  of  the  absorption 
coefficient  and  the  emissivity,  required  information  about 
the  spacing  of  the  individual  lines  within  the  group  and 
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their  relative  strength  is  supplemented  by  a  number  of 
models.  Two  extreme  models  are  Elsasser  model,  in  which 
equally  spaced  lines  of  equal  intensity  are  considered  and 
statistical  model,  in  which  the  spectral  lines  are  assumed 
to  have  random  spacing  and/or  intensity.  A  typical  spectral 
line  arrangement  for  these  models  is  shown  in  Figure  3. 


Figure  3.  Typical  Spectral  Line  Arrangements  for 

(a)  Elsasser  and  (b)  Statistical  Model  [1] 


In  statistical  model  it  is  assumed  that  the  spectral 
lines  are  not  equally  spaced  but,  rather,  are  randomly 
distributed  across  the  narrow  band.  A  statistical  analysis 
can  be  carried  out  in  two  different  ways.  In  the  uniform 
statistical  model  the  line  intensity  is  assumed  to  be 
constant,  leading  to 


W 

£  =  1  -  exp(— ) 
a 


where 


W 

d 


(2.5) 

-spectrally  averaged  emissivity 

-equivalent  line  width 
-spacing  between  spectral  lines 


In  general  statistical  model,  the  statistical  model 
has  also  been  evaluated  for  exponentially  decaying  line 
intensities,  resulting  in 


£n  =  1  -  CXP(_"T==)  =  1  “  exP(  ) 

i  T 

1  + -  1  +  — 


(2.6) 


where:  £  -spectrally  averaged  emissivity 


S  -line-integrated  absorption  coefficient 

d  -spacing  between  spectral  lines 

X  -optical  path  length 

bL  -line  half  width 


-line  overlap  parameter 

-narrow  band  optical  thickness 

The  results  from  statistical  model  is  summarized  in 
Figure  4  and  Table  1  in  [1]  as 
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Weak  line  Strong  line  No  overlap  All  regimes 


x«l 


line,  W  SX 


x  >  1 


0«1 


Elsasser  Model, 

S  -  const  1  - 1 
d  -  const 


Statistical  Model,  e  ^ 

S  =  const  1  -  e 
d  random 


S  exponential  1  -  e  7  1  -  expf- 


SXe  x[k M  +  /|(x)] 


2  TA[\-e~**M) 


er  T7y 


/  W 

l-exP-7 


d  random 


Definitions:  x  ~  fi  =  ir— ,  r  =  -X  -  2fix 
2vbL  d  d 
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Line  overlap  parameter  /3 


0.001 


0.0001  - 

- 1 - 1 _ _ lJ _ I _ 

0.001  0.01  0.1  1.0  10  100 

Line  strength  parameter  x 

Figure  4.  Regions  of  Validity  for  Elsasser  and 

Statistical  Band  Model  Regimes  [1] 


C.  RADIATIVE  TRANSPORT  EQUATION  (RTE) 

In  evaluating  the  radiative  energy  transport,  it  is 
required  to  know  radiative  properties,  such  as,  emissivity. 
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absorptivity,  and  reflectivity  in  the  case  of  surfaces,  and 
absorption  coefficient  and  scattering  coefficient  for  the 
participating  media  such  as  gases.  The  law  of  conservation 
of  energy  is  then  applied  to  determine  the  energy  field. 

The  general  relationship  that  governs  the  behavior  of 
radiative  heat  transfer  in  the  presence  of  an  absorbing, 

emitting  medium  is  developed  in  [1]  .  RTE  describes  the 
radiative  intensity  field  within  the  enclosure  as  a 

function  of  location,  fixed  by  location  vector  r, 

direction,  fixed  by  unit  direction  vector  ,  and  spectral 
variable,  wavenumber  .  The  net  radiative  heat  flux 

crossing  a  surface  element  is  obtained  by  adding  up  the 

contributions  of  radiative  energy  irradiating  the  surface 
from  all  possible  directions  and  for  all  possible 
wavenumbers.  Integrating  the  equation  of  transfer  over  all 
directions  and  wavenumbers  leads  to  a  conversation  of 
radiative  energy  statement  applied  to  an  infinitesimal 

volume . 

If  the  medium  is  participating,  in  that  case  any 
incident  beam  will  be  attenuated  by  absorption  and 
scattering  while  radiative  energy  travels  the  medium,  as 
shown  in  Figure  5. 

The  amount  of  absorption  is  directly  proportional  to 
the  magnitude  of  the  incident  energy  as  well  as  the 
distance  the  beam  travels  through  the  medium.  Thus, 

(dIv)abs  =  -K^nds  (2.7) 
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is  the  linear 


where  the  proportionality  constant  (7sn 

scattering  coefficient  for  scattering  from  the  pencil  of 
rays  under  consideration  into  all  other  directions. 

A  light  beam  traveling  through  a  participating  medium 
in  the  direction  of  s  loses  energy  by  absorption  and 
scattering  away  from  the  direction  of  travel.  But,  it  also 
gains  energy  by  emission  and  by  scattering  from  other 
directions  into  the  direction  s  . 

The  rate  of  emission  from  a  volume  element  is 
proportional  to  the  magnitude  of  the  volume.  The  emitted 
intensity  along  any  path  is  proportional  to  the  length  of 
the  path  and  the  local  energy  content  in  the  medium.  At 
thermodynamic  equilibrium,  the  intensity  everywhere  must  be 
equal  to  the  blackbody  intensity. 

(dln)em  =  KnIh11ds  (2.9) 

Augmentation  due  to  in-scattering  has  contributions 
from  all  directions,  and  it  is  calculated  by  integration 
over  all  solid  angles.  The  energy  flux  coming  from  all 
directions  and  scattered  into  the  direction  s  is  given  in 
[1]  as 

(dI„ ),«,  =dsJ^  j4„  ln  0"  (2.10) 

An  energy  balance  from  equations  (2.7),  (2.8),  (2.9), 

and  (2.10)  results  to  RTE  in  [1]  as 

ItJ(s+ds,s,t+ dt ) -In  (s,  s,  t )  =  KnIhn  (s,  t)ds  -  (s,  s,  t)ds  -  (s,  s,  t)ds 

a  ,  (2.11) 
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D. 


THE  METHOD  OF  DISCRETE  ORDINATES  (Sn-APPROXIMATION) 


The  Discrete  Ordinance  Method  (DOM)  is  a  tool  to 
transform  the  equation  of  transfer  into  a  set  of 
simultaneous  partial  differential  equations.  DOM  is  based 
on  a  discrete  representation  of  the  directional  variation 
of  the  radiative  intensity.  A  solution  to  the  transport 
problem  is  found  by  solving  the  equation  of  transfer  for  a 
set  of  discrete  directions  spanning  the  total  solid  angle 
range  of  4n  .  The  discrete  ordinates  method  is  simply  a 
finite  differencing  of  the  directional  dependence  of  the 
equation  of  transfer.  Integrals  over  solid  angle  are 
approximated  by  numerical  quadrature. 


The  general  equation  of  transfer  for  an  absorbing, 
emitting,  and  anisotropically  scattering  medium  is 


dl 

ds 


=  s  VI(r,s)  =  K(r)Ib{r)~ jB(r) I (r,s)  + 


4/r 


J  I(r,  s',  s)d£l' 

4  71 


(2.7) 


Equation  (2.7)  is  valid  for  a  gray  medium  or,  on  a 
spectral  basis,  for  a  nongray  medium,  and  is  subject  to  the 
boundary  condition 

I (rw ,  s)  =  £(rw )I b  (r  )  +  f  .  I(rw ,  i)|n  •  s\lCl'  (2.8) 

Jns<  0  1  1 

where  the  enclosure  is  considered  as  opaque,  diffusely 
emitting  and  reflecting  wall. 

In  the  discrete  ordinates  method,  equation  (2.7)  is 
solved  for  a  set  of  n  different  directions  sjf  i  =  l,2,0,n, 

and  the  integrals  over  direction  are  replaced  by  numerical 
quadratures,  that  is. 
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(2.9) 


£  f(s)dCl  = 

n  i=i 

where,  wt  is  the  quadrature  weights  associated  with  the 
directions  S;  . 

And,  equation  (2.7)  becomes, 

(7  ( r)  " 

■  V/(r,  s, )  =  K(r)Ib  ( r )  -  /3{r)I(r,  i, )  +  — j - ^  Wjl(r,  s  )4>(r,  i, ,  s  . ),  (2.10) 

4^  Jl 

where,  j=l,2,0,n,  subject  to  boundary  conditions 

I(rw,si)  =  £(rw)Ib(rw)  +  ^^-  (2.11) 

71  >’Sj<  o 

After  determining  the  intensities,  the  radiative  heat 
flux  inside  the  medium  or  at  a  surface  may  be  found  as 

n 

q(r )  =  [  I(r,  s)sdQ.  =  Y  wjIi  (r)st  (2.12) 

J47T  '  4 

i= 1 

In  Table  3.2,  the  direction  cosines  of  Sj  and 
corresponding  weights  for  S4 -approximation  is  shown.  The 
direction  cosines  of  Sj  are 

^  =  (s,  ■  i  )i  +  ( Si  ■  j)j  +  (s,  ■  k)k  =  %'i  +  rjJ  +  n ,i  (2.13) 

Only  positive  direction  cosines  are  given  in  Table 
3.2,  covering  one  eight  of  the  total  range  of  solid  angles 
4  .  To  cover  the  entire  4  any  or  all  of  the  values  of^,^, 
and  |i;  may  be  positive  or  negative. 

Only  positive  direction  cosines  are  given  in  Table  2, 
covering  one  eight  of  the  total  range  of  solid  angles  4  . 
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To 

and 


cover  the  entire  4  any  or  all  of  the  values  of^;, 
(.1;  may  be  positive  or  negative. 


Order  of 
Approximation 

Ordinates 

Weights 

5 

n 

(i 

w 

S4 

0.2958759 

0.2958759 

0 . 9082483 

0.5235987 

0.2958759 

0.9082483 

0.2958759 

0.5235987 

0 . 9082483 

0.2958759 

0.2958759 

0.5235987 

Table  2.  Discrete  Ordinates  for  the  S4 


Approximation  [1] 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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IV.  NUMERICAL  TOOLS 

A.  NUMERICAL  SOLUTION  STRATEGY 

There  are  mainly  three  numerical  tools  used  in  this 
study . 

First  tool  is  TEP,  which  is  used  to  calculate  the 
thermodynamic  chemical  equilibrium  properties  of  gas 
mixtures  resulting  from  the  combustion  of  various  fuels  and 
oxidants . 

The  second  one  is  RADCAL,  which  is  a  numerical  fire¬ 
modeling  program  which  provides  detailed  gas  spectral 
radiation  data. 

The  third  one  is  CFDACE,  which  is  a  set  of  computer 
programs  for  multi-physics  computational  analysis. 

The  numerical  solution  strategy  flowchart  is  shown  in 
Figure  6. 


Figure  6.  Numerical  Solution  Strategy  Flowchart 

B .  TEP  (THERMO-CHEMICAL  EQUILIBRIUM  PROGRAM) 

TEP  is  used  to  calculate  the  thermodynamic  chemical 

equilibrium  properties  of  gas  mixtures  for  a  variety  of 

different  applications.  Both  thermodynamic  data  and 

ingredient  libraries  are  included  with  TEP.  The 
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applications  that  are  included  in  the  TEP  program  are 
Rocket  Application,  Gas  Properties,  Shock  Application, 
Scramjet  Application,  Isentropic  Path,  and  Detonation 
Application . 

The  application  that  is  used  is  rocket  problem.  The 
rocket  problem  calculates  the  "theoretical",  or  maximum 
attainable  performance,  of  a  rocket  thrust  chamber.  The 
calculations  are  based  on  the  following  assumptions: 

•  One-dimensional  form  of  the  continuity,  energy,  and 
momentum  equations 

•  Zero  velocity  (stagnation)  in  the  combustion  chamber 

•  Complete  combustion 

•  Adiabatic  combustion 

•  Isentropic  expansion 

•  Homogeneous  mixing 

•  Ideal  gas  law 

•  Zero  temperature  and  velocity  lags  between  condensed 
and  gaseous  species 

For  equilibrium  performance,  composition  is  assumed  to 
reach  equilibrium  immediately  during  expansion.  For  frozen 
performance,  composition  is  assumed  to  remain  fixed  at  the 
initial  (stagnation)  combustion  composition  during 
expansion . 

Rocket  performance  is  obtained  by  first  determining 
stagnation  combustion  properties  in  the  rocket  chamber 
based  on  input  of  both  enthalpy  and  pressure  (H,  P)  or 
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temperature  and  pressure  (T,  P) .  Next,  conditions  at  the 
nozzle  throat  are  determined.  The  procedure  used  here  is 
to  find  the  pressure  at  which  the  flow  velocity  is  equal  to 
the  sound  speed.  Iteration  is  required  using  entropy, 
pressure  (S,  P)  point  solutions,  where  the  entropy  is  held 
fixed  at  the  value  calculated  previously  in  the  chamber. 
Once  the  throat  conditions  have  been  determined,  solutions 
at  input  area  ratios  are  obtained  by  iteration  using  (S,  P) 
point  solutions.  For  a  given  area  ratio,  A/A*,  the  entropy 
is  held  constant,  and  pressure  is  iterated  until  the 
product  of  density  and  velocity  satisfies  the  continuity 
relation;  i.e.,  a  value  equal  to  the  throat  density  x 
throat  velocity/  (A/A*). 

There  are  two  main  inputs  to  the  rocket  application, 
thermochemistry  and  rocket  operating  conditions.  In 
thermochemistry  chemical  reactants,  thermochemical  species 
data,  and  reactants  mixture  specification  is  defined.  Under 
thermochemistry  chemical  reactants  option,  the  ingredients 
that  will  be  used  as  oxidizer  and  fuel  are  selected.  And 
oxidizer  to  fuel  ratio  is  entered.  Under  thermochemical 
species  data  option,  the  species  that  will  be  used  in  the 
application  is  selected  between  CHNO,  CHO,  HO,  SRM,  Master 
Thermo,  and  Air/Hydrogen  file.  The  properties  of  these 
files  are  as  follows: 

CHNO  Thermodynamic  Data  File  contains  29  species 
screened  for  use  in  air/hydrocarbon  combustion  problems. 
This  set  of  species  is  more  than  sufficient  for  most  CHNO 
problems.  The  study  was  done  by  screening  calculated 

results  obtained  using  the  master  files.  The  calculations 
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were  for  air/acetylene  and  air/methane  for  a  broad  range  of 
pressure,  temperature,  and  equivalence  ratio. 

CHO  Thermodynamic  Data  File  contains  20  species 
screened  for  use  in  oxygen/hydrocarbon  combustion  problems. 
It  has  been  obtained  from  the  CHNO  thermo  data  file  by 
omitting  species  containing  N. 

HO  Thermodynamic  Data  File  contains  8  species  screened 
for  use  in  oxygen/hydrogen  combustion  problems.  It  has 
been  obtained  from  the  CHO  thermo  data  file  by  omitting 
species  containing  C. 

Solid  Propellant  Rocket  Motor  File  contains  the 
principal  species  that  make  up  the  exhaust  of  aluminized 
solid  propellant  rocket  motors. 

Master  Thermodynamic  Data  File  contains  data  for 
approximately  1000  chemical  species.  It  can  be  edited  to 
produce  short  files  containing  only  those  species  important 
to  specific  problems.  Several  files  of  this  type  are  also 
provided.  Principal  species  of  the  following  elements  are 


contained  in 

this  master 

f  ile : 

Al, 

Ar , 

B, 

Ba, 

Be, 

Br , 

c. 

Ca,  Cl,  Cr, 

Cs,  Cu,  D,  F, 

Fe,  H, 

He, 

Hg 

,  I,  K, 

Li, 

Mg, 

N, 

Na,  Nb,  Ne, 

Ni,  0,  P,  Pb, 

S,  Si, 

Sr, 

T  a , 

Ti, 

V, 

Xe, 

Zn, 

and 

Zr.  The  electron  and  positive  or  negative  ion  forms  for  a 
limited  number  of  species  are  also  provided.  Please  note 
that  the  master  file  is  intended  as  a  library,  and  not  for 
general  problem  solving  use.  The  master  file  will  often 
yield  too  many  species  for  problem  solving,  especially  when 
a  reactant  contains  the  element  carbon. 
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Air/Hydrogen  File  is  intended  for  use  in  Air/Hydrogen 
combustion  problems,  such  as  occur  in  flame  propagation  and 
shock  tube  studies. 

Under  rocket  operating  conditions  option,  chamber 
pressure  in  atmosphere,  chamber  temperature  in  Kelvin, 
subsonic  area,  supersonic  area,  and  pressure  ratios  are 
entered . 


Table  3  shows  the  input  and  output  data  for  TEP 
program . 


INPUT  DATA 

OUTPUT  DATA 

Gas  temperature,  gas 

Chamber  and  exit 

pressure,  oxidizer  to  fuel 

pressure,  participating  gas 

ratio,  subsonic  and 

temperature,  density. 

supersonic  area  ratios 

molecular  weight,  sonic 

velocity,  prandtl  number, 

mixture  species!  mole 

fractions,  mass  fractions. 

Table  3.  Input  and  Output  Data  for  TEP  Program 


C.  RADCAL 

RADCAL  is  a  narrow  band  numerical  model  developed  at 

NIST  for  radiation  calculations  in  a  combustion 

environment.  RADCAL  is  a  numerical  program,  which  predicts 

the  radiant  intensity  leaving  a  nonisothermal  volume 

containing  nonuniform  levels  of  carbon  dioxide,  water 

vapor,  methane,  carbon  monoxide,  nitrogen,  oxygen,  and 

soot.  The  absorption  coefficient  of  the  combined  gases  is 
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calculated  from  a  narrow-band  model  and  a  combination  of 


tabulated  spectral  properties  and  theoretical 
approximations  to  the  vibrational-rotational  molecular 
bands . 

Table  4  summarizes  the  species  and  their  molecular 
bands  currently  in  RADCAL. 


Species 

Band 

C02 

2.0  pm 

co2 

2.7  pm 

co2 

4.3  pm 

co2 

10.0  pm 

co2 

15.0  pm 

H20 

1.38  pm 

H20 

1.88  pm 

H20 

2.7  pm 

H20 

6.3  pm 

h2o 

20-200  pm 

CO 

4.6  pm 

ch4 

2 . 4  pm 

ch4 

3.3  pm 

ch4 

7.7  pm 

soot 

0.4-2000  pm 

Table  4.  Molecular  Bands  Included  in  RADCAL  [3] 

The  numerical  code  consists  of  a  short  main  program 
which  reads  the  temperature  and  concentration  information 
from  the  data  file,  RC  DAT.  SUBROUTINE  RADCAL  is  called 
from  the  main  program  to  perform  all  the  calculations. 
RCOUT  DAT  is  the  file  that  contains  the  outputs.  Figure  7 
is  a  diagram  of  the  structure  of  the  program. 

RADCAL  relies  upon  the  four  subroutines  to  compute  the 
narrow-band  parameters  for  the  carbon  dioxide 
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(SUBROUTINE  C02  )  ,  water  vapor  ( SUBROUT INE  H20  )  ,  carbon 
monoxide  (SUBROUTINE  CO )  ,  and  methane  (SUBROUTINE  FUEL).  Two 
large  block  data  files  (BLOCK  DATA  BD1  and  BLOCK  DATA  BD2 ) 
contain  the  absorption  coefficient  of  water  vapor  as  a 
function  of  temperature  and  wave  number.  A  third  data  file, 
BLOCK  DATA  BD3,  contains  similar  information  for  the  15  m 
band  C02  ,  and  the  3.3  and  7.4  m  bands  of  CHA. 

When  the  path  contains  only  one  element  three 
different  absorption  coefficients  are  calculated  the 
Planck-mean,  the  incident-mean,  and  the  effective 
absorption  coefficient. 

In  the  first  line  of  the  input  data  file,  RC  DAT,  the 
number  of  elements  into  which  the  path  is  divided,  NPT,  is 
written.  The  second  line  lists  the  size  of  the  first 
element  in  meters,  DD  (1),  its  temperature  in  Kelvin,  T 
(1),  and  the  partial  pressures  in  kilopascals  of  carbon 
dioxide,  P  (1,1),  water  vapor,  P  (1,2),  methane,  P  (1,3), 
carbon  monoxide,  P  (1,4),  oxygen,  P  (1,5), and  nitrogen,  P 
(1,6)  .  The  last  entry  on  the  second  line  is  the  volume 
fraction  of  soot  in  the  first  element,  W  (1)  .  The  third 
line  contains  the  size,  temperature,  species  partial 
pressure,  and  soot  volume  fraction  for  the  second  element. 

Similar  data  is  entered  for  the  remainder  of  the  NPT 
elements.  Following  this  information  is  a  line  containing 
the  wall  temperature  in  Kelvin,  TWALL,  and  the  minimum  and 
maximum  wavenumbers  incm"1,  OMMIN  and  OMMAX.  The  program 
terminates  when  it  reads  zero  for  the  value  of  NPT. 
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Figure  7.  Structure  of  Radiation  Calculation 

Program  RADCAL  [3] 


The  results  of  the  calculations  are  found  in  the 
output  file,  RCOUT  DAT.  The  input  conditions  are  summarized 
in  tabular  form,  followed  by  the  total  directional  radiated 
energy  flux  emanating  outward  from  element  one,  Q.  The 
spectral  intensity,  QW  (K)  and  the  transmittance,  TTAU  (K)  , 
are  listed  for  each  wavelength,  AMBDA  (K) .  The  program  also 
calculates  the  effective-mean  absorption  coefficient, 
AMEAN,  the  Planck-mean  absorption  coefficient,  APO,  and  the 
wall-incident-mean  absorption  coefficient,  AIWALL. 
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Table  5  shows  the  input  and  output  data  for  RADCAL 


program . 


INPUT  DATA 

OUTPUT  DATA 

The 

number  of 

elements , 

Total 

directional 

size  i 

of  the 

elements , 

radiated  energy 

flux,  the 

temperature  of  the 

elements , 

spectral  intensity  and  the 

partial 

pressures 

of  carbon 

transmittance 

for  each 

dioxide , 

water 

vapor. 

wavelength,  the 

effective- 

methane. 

carbon 

monoxide. 

mean  absorption 

coefficient  , 

oxygen. 

and  nitrogen,  volume 

the  Planck-mean 

absorption 

fraction 

of  soot,  wall 

coefficient,  and 

the  wall- 

temperature,  minimum  and 

incident -me an 

absorption 

maximum 

wavenumbers 

• 

coefficient . 

Table  5.  Input  and  Output  Data  for  RADCAL 


Program 


D.  CFD-ACE+ 

CFD-ACE+  is  a  set  of  computer  programs  for  multi¬ 
physics  computational  analysis  released  by  CFD  Research 
Corporation.  The  programs  provide  an  integrated  geometry 
and  grid  generation  module,  a  graphical  user  interface  for 
preparation  of  the  model,  a  computational  solver  for 
performing  the  simulation,  and  an  interactive  visualization 
program  for  examination  and  analysis  of  simulation  results. 
CFD-ACE+  package  include  the  following  applications: 
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•  CFD-GEOM  Interactive  3D  Geometry  Modeling  and 
Mesh  Generation  (Structured,  Unstructured,  and  Mixed- 
Element  Meshes) 

•  CFD-GUI  Solver  Setup  Interface 

•  CFD-ACE  (U)  General  Purpose  Multi-Disciplinary 
Physics  Solver 

•  CFD-VIEW  Interactive  3D  Graphics,  Animation  & 
Flow  Visualization  Software 


Figure  8  shows  the  applications  that  are  included  in 
CFD-ACE  +  package. 


CAD 


Figure  8.  CFD  ACE+  Package  [9] 

1 .  CFD-GEOM 

CFD-GEOM  is  an  interactive  CAD  type  geometry  creation 
and  fast  grid  generation  (structured,  unstructured  and 
hybrid  grids)  program  for  Computational  Fluid  Dynamics 
(CFD) .  The  geometry  module  allows  the  user  to  interactively 
define  the  geometry  or  modify  the  geometry  read  through  the 
Plot3D  data  format. 

It  serves  as  the  front-end  of  CFDRC ' s  solver  packages 


and  can  be  used  with  any  other  CFD/FEA/CAE  program.  With 
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interfaces  to  read  geometry  data  created  in  almost  any 
major  commercial  CAD  program,  it  integrates  CFDRC ' s  CFD 
packages  into  the  existing  CAE  environment. 

CFD-GEOM  uses  efficient  data  structures  providing 
flexible  inter-domain  connectivity,  part  merging,  and  easy 
manipulation  of  all  entities  (geometry,  topology, 
structured,  unstructured  meshes)  with  rapid  database 
updating.  Entity  manipulations  include  translation, 
rotation,  scaling,  copying,  (un) deleting,  blanking, 
labeling,  and  model  merging/extracting.  Topological 
entities  include: 

•  Edges  -  composed  of  single  or  multiple  curves  and/or 
lines . 

•  Faces  -  comprising  single  or  multiple  edge  sets  or 
single  surfaces,  possibly  projected  to  surfaces. 

•  Blocks  -  consisting  of  single  or  multiple  volumes  or 
face  sets.  Arbitrary  orientation  between  adjacent 
blocks  is  allowed,  or  a  consistent  multi-block 
orientation  can  be  requested  automatically  (hence  no 
need  to  pre-plan  individual  IJK  node  directions) . 

•  Composition  of  edges,  faces  and  blocks. 

•  Edge-linking  -  any  number  of  edges  can  be  linked  to  a 
smaller  number  of  master  edges  for  mesh  density  and/or 
distribution,  requiring  minimum  effort  to  modify  mesh 
characteristics . 

•  Boundary  Conditions-  the  location  and  type  of  boundary 
condition  for  the  flow  solvers  can  be  prescribed  in 
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CFD-GEOM,  providing  automatic  updating  during  mesh 

refinement  or  geometry  change. 

2 .  CFD-ACE  (U) 

CFD-ACE  is  an  unstructured,  polyhedral  cell  flow 
solver.  It  is  also  integrated  with  a  wide  variety  of  other 
physics  modules  making  it  the  core  of  a  multi-disciplinary 
analysis  environment.  CFD-ACE  employs  a  cell-centered 
control  volume  solution  approach.  This  approach  implies 
that  the  discrete  equations  are  formulated  by  evaluating 
and  integrating  the  fluxes  across  the  faces  that  surround 
each  control  volume. 

CFD-ACE  provides  an  integrated  geometry  and  grid 
generation  module. 

3 .  CFD-VIEW 

CFD-VIEW  is  a  3-D  post  processor  included  in  CFD-ACE+, 
which  reduces  the  large  volume  of  data  generated  by  CFD-ACE 
to  useful  information.  CFD-VIEW  is  an  interactive  graphics 
program  for  post-processing  numerical  results  from  CFD  and 


other  analysis 

software . 

It 

provides 

an  easy- 

-to-use 

and 

interactive  environment. 

with  many 

graphics 

tools 

to 

visualize  the 

flow  physics. 

as  well 

as  the 

ability 

to 

extract  data  relevant  to  engineering  design. 
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V. 


CODE  CALIBRATION  AND  BENCHMARKING 


A.  NUMERICAL  TOOL  CALIBRATION 

Geometry  of  the  model  is  created  by  the  pre-processor 
CFD-GEOM.  Creating  the  structured  grid  is  done  by  grid 
generation  module  in  a  bottom-up  approach  in  labeling 
geometry  elements  as  edges,  faces,  and  blocks.  In  order  to 
achieve  the  best  results,  the  calibration  is  done  when  the 
participating  medium  is  not  present.  This  helps  to  compare 
the  results  with  the  hand  calculations. 

The  following  parameters  are  considered  when 
calibrating  the  model: 

•  The  edge  parameters 

•  The  number  of  iterations 


•  The  error  in  the  wall  convection 


The  model  is  created  by  25,  50,  100,  and  200  points 
consecutively.  The  edge  parameters  dictate  how  many  grid 
points  are  placed  on  the  edge  and  how  these  grid  points  are 
distributed.  The  power  law  is  used  as  a  grid  point 
distribution  method.  The  power  law  distribution  is  defined 
by 


u(n)  = 


f 

v 


n  —  1 
npts  - 1 


y 

) 


(5.1) 


where:  u(ri)  -location  of  grid  point  on  the 

interval  [ 0  0 1 ] 
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n  -grid  points  index  on  the  interval 

[  lOnpt s ] 

npts  -the  number  of  grid  points 


x  -user-specified  power 

The  power  law  is  usually  close  to  one  which 
corresponds  to  a  uniform  distribution.  For  non-unity 
powers,  the  forward,  backward,  and  symmetry  options 
indicate  how  the  power  law  grid  distribution  is  applied. 


The  geometry  is  created  with  25,  50,  100  and  200 

points,  while  the  power  law  grid  scaling  scheme  is  varied 
between  1.3,  1.5  and  2.  Convergence  criteria  were  varied 

from  le-04  to  le-10.  Figure  9  discusses  the  results  for 
convergence  criteria  of  le-08. 
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Figure  9.  Total  Heat  Flux  Values  for  Various 

Power  Law  Grid  Scaling  Schemes 

As  seen  in  the  figure,  the  heat  flux  value  for  100 
points  with  power  law  of  2  converges  more  rapidly  than  the 
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others  and  a  100  point  grid  was  found  to  be  a  good  trade¬ 
off  between  accuracy  and  computational  time. 

B.  TRASPARENT  MEDIUM  CASE 

Numerical  results  are  compared  with  the  analysis  of 
radiative  heat  transfer  within  enclosures  without 
participating  medium.  Figure  10  defines  the  surfaces  and 
the  emissivities  corresponding  to  these  surfaces. 


p  T 


Figure  10.  Two  Dimensional  Cylinder 

Radiative  exchange  between  black  surfaces  is 
calculated  from  the  following  expression  in  [1]  as 

9t=ftF,-j(Ebi-Eli)  (5.2) 

i=i 

where : 

N  -number  of  surfaces 

Fj  j  -view  factor  from  surface  i  to  j 

Ehi  -blackbody  emissive  power  of  surface  i 

The  equation  5.2  becomes 
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‘h  ~  ^1-1  EbX ) +  F\-i  (Em  Eb2 )  +  ^i_3  (Ehl  Eb3 ) 


(5.3) 


where  the  view  factors  can  be  calculated  from  view  factor 
catalogue.  [2] 


Table  6  shows  the  comparison  of  the  hand  calculations 
and  CFD-A  software  results  for  wall  temperature  of  400  F 
and  gas  temperature  of  3500  K  between  black  surfaces. 


Hand 

Calculation 

Numerical 

Results 

Relative  Error 

(%) 

Radiation  to  the 

Wall  (watts/rad) 

-2882.46 

- 

3.02 

2795 

Table  6.  Relative  Error  in  Radiative  Exchange 

Between  Black  Surfaces 


Similarly,  radiation  exchange  between  gray  and  diffuse 
surfaces  can  be  calculated  from  the  following  expression 
[1] 


c1l 


N  1 

K--1^ 


2 ]Fi-j(Ebi  ~Ebj) 


7=1 


(5.4) 


The  equation  (5.4)  becomes 

—  -  “  ^Wx-x^x  =  Fx-x  (Ebl  ~Ebl)  +  Fx_2 (Ebl  ~Eb2)  +  Fx_3  (Ebl  -Eh3)(  5.5) 

£x  £i 
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Table  7  shows  the  comparison  of  the  hand  calculations 
and  CFD-ACE  software  results  for  wall  temperature  of  400  F 
and  gas  temperature  of  3500  K  between  gray  surfaces. 


The  relative  error  is  due  to  the  fact  that  in  hand 
calculations  the  outlet  is  assumed  to  be  isothermal,  which 
is  not  the  case  for  real  cases. 


Hand 

Calculation 

Numerical 

Results 

Relative 

Error  (%) 

Radiation 

-2512.24 

-2369 

5 .701 

to  the  Wall 

(watts/ rad) 

Table  7.  Relative  Error  in  Radiative  Exchange 


Between  Gray  Surfaces 


C.  GRAY  GAS  CASE 

A  literature  search  shows  that,  although  there  is  an 
extensive  search  for  the  gray  gas  approximation,  the 
pressure  and  temperature  values  are  limited  to  the 
temperatures  and  pressure  values  lower  than  that  of 
interest  in  this  study.  [2, 6, 7, 8]  .  Next,  CFD-ACE  results  are 
compared  with  hand  calculations  for  different  optical 
thickness  values  for  parallel  plate  and  round  tube  cases 
with  participating  medium.  Parallel  plate  and  round  tube 
cases  are  considered  consecutively. 
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1. 


Parallel  Plate 


In  parallel  plate  case,  a  gray,  nonscattering  medium 
is  contained  between  two  parallel  gray  plates.  The  medium 
is  isothermal  at  temperature  Tg  ,  with  constant  absorption 

coefficient  k  .  The  two  plates  are  isothermal  at  temperature 
Tw ,  have  the  same  gray-diffuse  emissivity  ,  and  are  spaced 
a  distance  D  apart  as  shown  in  Figure  11. 


Figure  11.  Coordinate  System  and  Dimensions  for 

the  Parallel  Plate 

Heat  flux  at  for  the  specified  conditions  is  given  in 
[1]  as 


q(T)  =  cr(Tn4 


T  4 \  2 [E3(t)-E3(tl-t)\ 

l  +  (--l)[l-2  E3(tl)\ 
£ 


(5.6) 


At  the  wall  x  is  equal  to  zero.  So,  the  equation 
becomes 


q(T  =  0)  =  a(Tw-T;) 


4  m  4\  2[£’3(0)  E3(tl )] 


l  +  (--l)[l-2  E3(tl)\ 
£ 


(5.7) 
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and 


q(r  =  0)  =  g (Tj  -  Tg 4) - [1  E3(Ti)] -  (5.8) 

1  +  (— 1)[1-  2E3(tl)] 

£ 

Equation  5.8  is  valid  for  all  the  optical  thickness 
values  other  than  zero  and  infinity.  For  optically  thick 
and  optically  thin  cases,  two  different  solutions  are 
necessary . 

In  optically  thick  case  xL  goes  to  infinity  and 
Equation  (5.8)  becomes 

q  =  £G(Tw4-T g4)  (5.9) 

In  the  optically  thin  situation  the  following  result 
for  total  heat  loss  is  used  for  any  isothermal  volume 
without  self  absorption  in  [1] 

Q  =  4m2oTg4V  (5.10) 


Multiplying  and  dividing  this  equation  with  2R  and 
writing  the  heat  flux  gives 

o  4rrn2oT4V/2R  2rfn2oT4 
q  -  As  ~  (2 m  +  2 Di  +  2 HD)  ~  1  D_  D 

H  i 


where,  V  is  total  volume  and  As  is  the  surface  area. 

Assumption  of  l  is  greater  than  D  and  H  is  a  very  good 
one  and  equation  (5.11)  becomes 


q  = 


4TLn2cfiyi2R  2  TLn2oTg 


Q__ _ _ _ 

As  (2  Hi  +  2  Di  +  2  HD)  }  D_ 

H 


(5.12) 
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Since  the  case  is  optically  thick,  taking  two  values 
for  tl  is  good  enough,  such  as  zero  and  one.  The  other 
unknown  is  D/H  ratio  and  for  this  values,  zero  and  one  and 
a  half  is  chosen.  Zero  corresponds  to  a  very  long  parallel 
plate,  which  makes  H  equal  to  infinity,  but  it  causes  to 
deteriorate  the  true  geometry.  On  the  other  hand,  one  and  a 
half  corresponds  to  the  real  geometry.  Figure  12  and  13 
compare  the  theoretical  radiation  to  the  wall  values  with 
CFD-ACE  results  for  various  optical  thicknesses,  for 
optically  thick  and  optically  thin  case  in  a  parallel 
plate . 


Figure  12.  Comparison  of  Theoretical  and  Numerical 

Results  for  Parallel  Plate 
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Figure  13. 


Comparison  of  Theoretical  and  Numerical 
Results  for  Parallel  Plate  (magnified) 


2.  Round  Tube 


In  round  tube  case,  a  gray,  nonscattering  medium  is 
contained  between  the  walls  of  the  cylinder.  The  medium  is 
isothermal  at  temperature  Tg  ,  with  constant  absorption 

coefficient  k.  The  wall  of  the  cylinder  is  isothermal  at 
temperature  Tw ,  have  the  gray-diffuse  emissivity  ,  and  the 
radius  is  R  as  shown  in  Figure  14. 

The  nondimensional  heat  loss  from  a  gray, 
nonscattering,  isothermal  cylinder  is  given  in  [1]  as 


w  = 


4(?r) 


nz<trg-Jn 


(5 . 13) 
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Since  the  walls  are  not  black,  but  are  gray,  diffuse 
emitters  and  reflectors,  the  heat  flux  to  the  wall  is  no 


longer  7ilb ,  but  must  be  replaced  by  the  radiosity. 


q(T„ )  =  Y(n-<JT;  -  JJ  =  (. J „  -E„) J- 

1  -  £ 


(5.14) 


Since  the  walls  are  not  black,  but  are  gray,  diffuse 
emitters  and  reflectors,  the  heat  flux  to  the  wall  is  no 
longer  7ilb ,  but  must  be  replaced  by  the  radiosity. 


q(TR)  =  ¥(,rcjr;  -JJ  =  C Jw-EbJ  £ 


l-£„ 


So,  the  radiosity  becomes 


J , 


(^X4)  +  (-^orw4) 


1  ~£ 


+  ¥ 


And  the  heat  flux  is 


q(rj  =  ay/ 


(n%4)- 


(¥n2T  4)  +  (-£^Tw4) 

1  - 


1  -£ 


+  ¥ 


(5.14) 


(5 . 15) 


(5.16) 


Equation  5.16  is  valid  for  all  the  optical  thickness 
values  other  than  zero  and  infinity.  For  optically  thick 
and  optically  thin  cases,  two  different  solutions  are 
necessary . 

In  optically  thick  case  v|/  goes  to  one  according  to 
Table  8  and  Equation  (5.16)  becomes 
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r)  —  d" 


(5.17) 


(n2r/)  +  (-^rw4) 

(«V) - — 


In  the  optically  thin  situation  the  following  result 
for  total  heat  loss  is  used  for  any  isothermal  volume 
without  self  absorption  in  [1] 

Q  =  4m2oTg4V  (5.18) 


Multiplying  and  dividing  this  equation  with  R  and 
writing  the  heat  flux  gives 


q  = 


Q 

A 


4 lRH  ^  VIR  =  2TRn2&r  \—[— 

(ttR2  +ttR2  +  2tiRH)  s  Rl  +  1 

H 


(5.11) 


where,  V  is  total  volume  of  the  cylinder  and  As  is  the 
surface  area. 


Since  the  case  is  optically  thick,  taking  two  values 
for  tr  is  good  enough,  such  as  zero  and  one.  The  other 
unknown  is  R/H  ratio  and  for  these  values,  zero  and  one 
quarter  is  chosen.  Zero  corresponds  to  a  very  long 
cylinder,  which  makes  H  equal  to  infinity,  but  it  causes  to 
deteriorate  the  true  geometry.  On  the  other  hand,  one 
quarter  corresponds  to  the  real  geometry.  Figure  15  and  16 
compares  the  theoretical  radiation  to  the  wall  values  with 
CFD-ACE  results  for  various  optical  thicknesses,  for 
optically  thick  and  optically  thin  case  in  a  round  tube. 
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radiation  to  the  wall  (W/rad) 


Figure  15.  Comparison  of  Theoretical  and  Numerical 

Results  for  Round  Tube 


optical  thickness 

Figure  16.  Comparison  of  Theoretical  and  Numerical 

Results  for  Round  Tube  (magnified) 
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It  is  obvious  that,  as  the  participating  gas  gets 
optically  thick  the  results  do  not  agree  since,  in  reality, 
the  gas  is  not  isothermal  across  the  chamber  cross-section. 

When  the  gas  is  optically  thick,  its  absorption 
coefficient  is  large,  so  the  typical  absorption  length 
scale  is  small.  In  other  words,  emission  from  the  gas  is 
reabsorbed  over  a  very  short  distance.  In  such  a  scenario, 
the  radiation  on  the  wall  comes  from  the  region  of  gas 
close  to  the  wall  where  the  gas  is  at  a  lower  temperature 
than  the  core  temperature  (due  to  the  boundary  layer  nature 
of  the  temperature  profile) ,  and  is  hence  a  lower  heat  flux 
value . 

An  important  point  is  to  investigate  if  these  hand 
calculations  contain  radiation  to  the  wall  according  to  the 
expressions  used  herein  or  not.  So,  in  order  to  compare 
these  results  with  each  other,  the  value  of  radiation  to 
the  wall  for  gray  walls,  transparent  medium  from  CFD-ACE, 
is  not  added  to  the  hand  calculation  results.  This  is  shown 
in  Appendix  A. 

D.  NON-GRAY  GAS  CASE 

A  literature  search  shows  that  only  a  limited  number 
of  non-gray  cases  have  been  treated  in  the  past,  and  have 
been  restricted  to  pure  gases  at  temperatures  and  pressures 
lower  than  that  of  interest  in  this  study  [4,5]  .  The 

radiation  heat  flux  in  CO  for  1000  °R  gas  temperature  at  1 
atmosphere  pressure  and  500  °R  wall  temperature  with  a 
parabolic  initial  velocity  distribution  is  considered.  The 
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parabolic  velocity  profile  at  the  inlet  is  created  by  using 
the  following  expression  and  shown  if  Figure  17. 


U  =  2Uave(l- 


(5.12) 


Figure 

Where, 


U 


17  . 

Uave  iS 

 Re*  ,4/ 

ave  p. 

p*D 


Parabolic  Velocity  Profile  at  the  Inlet 
found  from  the  following  expression. 

(5 . 13) 


By  using  the  fire  prediction  program  RADCAL,  the 
intensities  and  transmissivities  are  found  based  on  the 
wavelength  as  shown  below  in  Figure  18. 

The  numerical  result  and  published  result  are  compared 
in  Figure  19  and  compared  with  [4]  .  The  published  results 
and  numerical  results  agree  well. 
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Rq/kT  transmissivity 


8.  RADCAL  Output  for  Transmissivity  of 


Pure  CO 


1.E-04  1.E-03  1.E-02  1.E-01 


(x/R)/RePr 


igure  19. 


CFD  Output  for  the  Specified  Conditions 


VI.  RESULTS  AND  DISCUSSIONS 


The  model  problem  being  treated  in  this  study  contains 
several  parameters,  such  as  the  gas  optical  thickness,  the 
conduction  to  radiation  parameter,  the  emissivity  and 
temperature  of  the  wall,  the  gas  mixture  constituents  and 
its  temperature,  etc.  Of  these,  only  the  parameters  of 
primary  importance  were  varied  in  a  parametric  study  as 
described  below  so  as  to  determine  the  role  of  the  main 
contributing  factors  to  the  radiative  heat  loading  of  the 
duct  wall . 

The  results  presented  here  are  obtained  using  a 
100x100  grid.  Further  subdivision  of  the  enclosure  does  not 
yield  significantly  more  accurate  results  under  the 
specified  convergence  criteria.  Convergence  is  checked 
using  a  tolerance  of  le-10.  The  CPU  time  required  to 
generate  the  results  varied  from  22  minutes  to  107  minutes 
(per  run)  depending  on  the  complexity  of  the  case. 

A.  GRAY  GAS  CASE 

Under  the  gray  gas  approximation,  a  parametric  study 
was  conducted  by  varying  gas  temperature,  wall  temperature, 
wall  emissivity,  and  optical  thickness.  The  role  of  the  gas 
temperature  for  the  gray  gas  approximation  is  shown  in 
Figure  20  and  21.  Figure  20  shows  the  convective  and 
radiative  heat  fluxes  to  the  wall  versus  the  optical 
thickness  for  a  low  gas  to  wall  temperature  ratio  and  for  a 
high  gas  to  wall  temperature  ratio.  It  is  evident  that 
there  is  an  optimal  optical  thickness  at  which  the 


47 


radiative  flux  is  a  maximum,  which  can  be  seen  clearly  in 
Figure  20.  As  expected  the  radiative  flux  is  higher  on  the 
black  wall,  and  higher  gas  temperatures.  In  both  figures, 
it  is  shown  that,  the  radiative  flux  is  a  significant 
fraction  of  the  total  heat  flux  especially  in  the  large 
gas-to-wall  temp  ratio  case. 
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Figure  2  0 . 


Convective  and  Radiative  Heat  Flux  for 


1500K  Gas  800K  Wall  Temperature 


The  dependence  on  wall  emissivity  for  the  gray  gas 
approximation  is  shown  in  Figure  22  and  Figure  23.  These 
figures  show  the  relative  contribution  of  the  radiative 
heat  flux  versus  the  non-dimensional  length  along  the  duct 
for  low  gas  to  wall  temperature  and  for  high  gas  to  wall 
temperature  ratios,  respectively.  An  optical  thickness  of  3 


48 


is  used  which  is  close  to  the  optimal  value  as  shown  in  the 
previous  figures  for  maximum  radiative  flux.  It  is  clear 
that  the  relative  radiative  contribution  for  the  high  gas 
temperature  is  almost  twice  that  for  the  low  temperature 
case . 


0  2  4  6  8  10 

optical  thickness 

Figure  21.  Convective  and  Radiative  Heat  Flux  for 

1500K  Gas  800K  Wall  Temperature 

Figures  24  and  Figure  25  show  the  radiative  heat  flux 
versus  the  non-dimensional  length  along  the  duct  for 
different  wall  temperatures  at  3500  K  gas  temperature. 

It  is  clear  in  Figure  25  that,  the  absolute  heat  flux 
values  are  nearly  the  same,  although  in  Figure  24,  the 
relative  contribution  is  slightly  higher  for  the  hotter 
wall . 
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Figure  22.  Relative  Contribution  of  the  Radiative 

Heat  Flux  in  Different  Wall  Emissivities 
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Figure  23.  Relative  Contribution  of  the  Radiative 

Heat  Flux  in  Different  Wall  Emissivities 
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Figure  24.  Relative  Contribution  of  the  Radiative 

Heat  Flux  in  Different  Wall  Emissivities 
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Figure  25.  Relative  Contribution  of  the  Radiative 

Heat  Flux  in  Different  Wall  Emissivities 
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B . 


NONGRAY  GAS  CASE 


A  parametric  study  was  conducted  by  varying  gas 
pressure  in  addition  to  the  parameters  which  used  in  the 
gray  gas  approximation.  Figure  26  shows  the  typical 
constituents  in  the  gas  mixture  resulting  from  the  complete 
combustion  of  RP1  fuel  with  oxygen  as  oxidant  and  an  0/F 
ratio  of  3.0.  Figure  26  is  for  the  particular  gas  of  a  1500 
K  gas  mixture  at  a  pressure  of  25  atmospheres. 


Figure  26.  Typical  Mole  Fraction  of  A  Gas  Mixture 


These 

mole  fractions  are 

used  as 

input 

in  RAD CAL 

and 

to  obtain 

spectral  radiation 

data  as 

shown 

in  Figure 

27 

where  transmissivity  is  plotted  versus  wavelength.  The 
bandwidth  of  each  constituent  is  divided  into  wavenumbers. 

The  Figures  28  and  29  discuss  inlet  and  outlet 
emissivity  and  show  the  radiative  heat  flux  to  the  wall  for 
low  gas  to  wall  temperature  ratio,  and  for  high  gas  to  wall 
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temperature  ratio  at  a  chamber  pressure  of  35  atmosphere. 
The  blue  line  is  for  black  inlet  &  outlet  surfaces,  for 
which  the  radiative  heat  flux  is  a  maximum. 


Figure  27.  Typical  Spectral  Transmissivity 
Variation  of  the  Mixture 

The  purple  line  is  for  the  case  of  zero  inlet/outlet 
emissivity  and  shows  only  the  core  gas  contribution.  The 
0.5  emissivity  curve  is  an  intermediate  result  between  the 
previous  cases. 

Figures  30  and  31  discuss  the  role  of  wall  emissivity 
by  comparing  black  and  gray  walls  and  show  the  relative 
contribution  of  the  radiative  heat  flux  versus  the  non- 
dimensional  length  along  the  duct  for  low  gas  to  wall 
temperature  ratio  and  high  gas  to  wall  temperature  ratio  at 
a  chamber  pressure  of  35  atmosphere  pressures. 
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It  is  clear  that,  the  heat  flux  is  about  50%  higher  to 
black  walls,  which  appears  to  be  fairly  independent  of  the 
actual  gas  temperature. 


(x/R)/RePr 


Figure  28.  Heat  Flux  for  High  Gas-to-wall 

Temperature  Ratio 


1.E-06  1.E-05  1.E-04 

(x/R)/RePr 


Figure  29.  Heat  Flux  for  Low  Gas-to-wall 

Temperature  Ratio 
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1.E-06  1.E-05  1.E-04 

(x/R)/Re  Pr 

Figure  30.  Heat  Flux  for  High  Gas-to-wall 

Temperature  Ratio  for  Black  Wall 


1.E-06  1.E-05  1.E-04 

(x/R)/RePr 


Figure  31.  Heat  Flux  for  High  Gas-to-wall 

Temperature  Ratio  for  Gray  Wall  (emissivity=0 . 5 ) 

Figures  32  and  33  discuss  role  of  gas  pressure  at  3500 
K  gas  temperature  and  400  K  wall  temperature  for  the  case 
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of  a  black  wall  and  show  the  variation  of  heat  flux  versus 


non-dimensional  length  for  various  chamber  pressures.  Both 
the  absolute  and  relative  values  of  the  radiative  heat  flux 
do  not  change  appreciably  due  to  the  fact  that  the 
absorption  coefficients  in  the  most  active  bands  are 
already  around  their  optimal  value,  and  so  much  larger 
pressures  are  needed  to  alter  the  radiative  heat  flux 
significantly . 


1.E-06  1.E-05  1.E-04 

(x/R)/RePr 

Figure  32.  Relative  Contribution  of  Heat  Flux  in 

Black  Wall  for  Various  Gas  Pressures 


1.E-06  1.E-05  1.E-04 

(x/R)/RePr 

Figure  33.  Radiative  Heat  Flux  in  Black  wall  for 

Various  Gas  Pressures 
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VII .  CONCLUSIONS  AND  RECOMMENDATIONS 


A  numerical  model  has  been  developed  for  finding  the 
radiative  heat  load  on  the  walls  of  a  combustion  chamber. 
Such  a  model  is  useful  for  improving  and  augmenting  the 
predictive  capabilities  of  other  numerical  codes  that 
provide  only  convective  heat  loads. 

Although  some  of  the  trends  in  results  are  expected, 
this  study  provides  concrete  quantitative  estimates  of  the 
role  of  gas  temperature,  wall  temperature  and  emissivity, 
etc . 

It  was  found  that  the  radiative  fluxes  are  least 
sensitive  to  the  wall  temperature  and  chamber  pressure  in 
the  range  considered.  However,  the  effective  emissivities 
of  the  inlet  and  outlet  surfaces  were  found  to  be  crucial 
in  determining  the  magnitude  of  these  fluxes. 

Most  importantly,  it  was  found  that,  a  proper  non-gray 
spectral  treatment  of  the  gas  mixture  is  critical  for  an 
accurate  analysis.  A  gray  treatment  is  too  simplistic  and 
yields  unreliable  results. 

More  accurate  and  realistic  estimates  can  be  obtained 
by  expanding  the  study  to  include  the  following: 

•  The  absorption-emission  and  scattering  effects  of 
particulate  matter  such  as  soot 

•  The  role  of  the  liquid  fuel  that  is  injected  to  form  a 
film  on  the  chamber  walls 
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•  The  interaction  of  the  combustion  mechanics  with  the 


flow  and  heat  transfer  characteristics 

•  The  effect  of  the  resulting  non-uniform  distribution 
of  gaseous  species  in  the  chamber 
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APPENDIX  A.  THE  WALL  CONTRIBUTION  TO  HEAT  FLUX 


If  the  medium  is  not  at  radiative  equilibrium  (i.e. 
conduction  and/or  convection  is  not  negligible  and 
radiative  heat  transfer  is  not  dominant)  the  heat  flux  is 
found  by  using  the  definition  of  the  exponential  integrals 
in  [1]  as 

T  rL 

q(T)  =  2 JxE3  (t)  -  2  J2E3  (tl  -  r)  +  2oTm4  J  E2  (r  -  -  2 oTm4  J  E2  (t  -  r)dr 

0  T 

and  it  becomes 

q(T)  =  2JlE2(T)-2J2Ei(TL-T)  +  2oTm4[E2(T)-Ei(TL  ~  T)] 


Figure  34.  Optical  Coordinates  for  the  Model 

As  shown  in  Figure  34  optical  coordinate  at  the  upper 
wall  is  zero  and  optical  coordinate  at  the  lower  wall  is  xL . 

Setting  x  equal  to  zero  and  xL  in  the  previous  equation 
consecutively  gives 

<7(0)  =  2  J,E3  (0)  -  2  J2E2  (tl  )  +  2 oTm 4 [E 3  (0)  -  E,  (rL )]  =  W 4  -  J, ) 

l-^i 

q(rL)  =  2J \E3 (tl ) - 2 J2E, (0)  +  2&Tn; [E3 (r) - E 3 (0)]  =  ( oT2 4  -  J2 ) 

\-e2 
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As  seen  in  the  equation  above,  the  expression  found 
for  the  heat  flux  at  the  upper  wall  considers  the 
emissivity  and  temperatures  of  the  top  and  bottom  walls. 
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APPENDIX  B.  TYPICAL  CFD-ACE  OUTPUT 


-k  ■k  ■k  ★ 


★  ■k 

CCCCC  FFFFF 

DDDD 

AAA 

CCCCC 

EEEEE 

U 

u 

■k  ■k 

C 

F 

D  D 

A  A 

C 

E 

U 

u 

■k  ★ 

k  it 

c 

FFFF 

D  D  = 

=  AAAAA 

C 

EEEE 

U 

u 

■k  k 

k  ■k 

c 

F 

D  D 

A  A 

C 

E 

U 

u 

■k  k 

k  ■k 

CCCCC  F 

DDDD 

A  A 

CCCCC 

EEEEE 

UUU 

■k  k 

k  ■k 

■k  ★ 

k  it 

Version 

2002.  0. 

27 

■k  ■k 

Build  Date  : 

07/28/2002  22:3 

4  :  49 

k 

k  ■k 

Build  OS 

Windows_ 

NT 

■k  ★ 

k  ■k 

Build  OS 

Release  : 

1.3.12  (0 

.54/3/2 

) 

■k  k 

★  ■k 

Build  OS 

Version  : 

2002-07- 

06  02:1 

6 

■k  k 

k  ■k 

Build  Machine  : 

BELL2 

■k  k 

k  ■k 

■k  k 

k  ■k 

Copyright 

(c)  2000, 

CFD  Research  Corporat 

ion. 

■k  ★ 

★  ■k 

All  Rights  Reserved 

■k  ★ 

CFD-ACEU  Run  Platform  Information  : 


Run 

Date 

:  11/24/2002 

Run 

OS 

:  Windows 

Run 

OS  Release 

Run 

OS  Version 

Run 

Machine 

:  IT002099 
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Zone  No. 


1 


VC  Name .  :  NoName 

No.  of  Cells  :  9801 

Material  Type  :  Fluid 

Density  method  :  Constant  =  2.97E+00 

Viscosity  method  :  Constant_Dyn  =  9.52E-05 

Cond.  method:  Prandtl  =  6.5E-01 

Sp .  Heat  method:  Constant  =  1.94E+03 

Absorption  Coeff.  Set:  gas 
Emissivity  Set:  None 


RAD_SNORD  Statistics 


Parameters 


Radiation 
Gray  model 
CASN-DOM 

S4  scheme 
sn_hrange  scheme 
%  Relative  Error 
err_0m 
err_lm_mu 
err_lm_xi 
err_lm_et 
err_2m_mu 
err_2m_xi 
err_2m_et 


=  TRUE 

=  TRUE 

=  TRUE 

=  TRUE 
in  Moments 
-2 . 980232E-06 
-1 . 092752E-05 
-1 . 092752E-05 
-1 . 000000E+02 
-1 . 185405E-05 
-1 . 185405E-05 
-1 . 185405E-05 
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Total  Mass  Flow  Summary 


2 . 786166E-02 


2 . 786166E-02  -5.55112E-17 
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Name 


Key  Type 


X-axis 


Y-axis 


Z-axis 

1  NoName  28  Wall  0.000000E+00  0.000000E+00 

1 . 778739E+02 


Viscous  Moments 


Name 


Z-axis 


1  NoName 
-1 . 392094E-03 


Key  Type  X-axis  Y-axis 


28  Wall  0.000000E+00  0.000000E+00 


Boundary-by-Boundary  Heat  Transfer  Summary  (watts/rad) 


Name 

key 

Type 

COND . +CONV . 

W_SRC . +CVD 

RADIATION* 

Sum 

NoName 

37 

Inlet 

1 . 8899E+05 

0.0000E+00 

1 . 2769E+02 

1 . 8912E+05 

NoName 

36 

Outlet 

-1 . 7  693E+05 

0.0000E+00 

3 . 6525E+02  - 

1 . 7656E+05 

NoName 

28 

Wall 

-8 . 2173E+03 

0.0000E+00 

-4 . 3441E+03  - 

1 . 2561E+04 


Total  volume  source 
0 . 0000E+00 
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Total  Heat  Imbalance  3.8479E+03  0.0000E+00  -3.8511E+03  - 

3 . 2126E+00 


Total  wall  HEAT  source 
0 . 0000E+00 

*  Radiation  Summary  Convention  is  +ve  if  emmitting,  -ve  if  absorbing. 

*  Cond+Conv  Summary  Convention  is  +ve  if  flux  into  the  cell,  -ve  if 
flux  leavi  ng  the  cell. 


End  of  Iterative  Cycle 


Final  Time  Elapsed  Time=  6.552923E+02  Delta-time=  6.552923E+02 
Normal  Termination 
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APPENDIX  C.  TYPICAL  RADCAL  INPUT 


l 


0.2850 

1000.  7.5  15.6 

.  0 

.  0 

.  0 

77 . 0 

.  0 

0.  50. 

1 

10000  . 

1 . 000 

1000.  7.5  15.6  . 

.0  , 

.0  , 

.0 

77.0  , 

.  0 

0.  2000.  3000. 

0 

0.030 

1000.  7.5  15.6  . 

.0  , 

.0  , 

.  0 

77.0  , 

.  0 

1500  . 

50.  10000. 

1 

0.010 

1000.  7.5  15.6  . 

.0  , 

.0  . 

.  0 

77.0  . 

.  0 

1500  . 

50.  10000. 

1 

0.003 

1000.  7.5  15.6  . 

.0  , 

.0  , 

.  0 

77.0  , 

.  0 

1500  . 

50.  10000. 

1 

0.001 

1000.  7.5  15.6  . 

.0  , 

.0  , 

.  0 

77.0  , 

.  0 

1500  . 

50.  10000. 

1 

0.0003 

1000.  7.5  15.6 

.  0 

.  0 

.  0 

77 . 0 

.0 

1500  . 

50.  10000. 

1 

0.0001 

1000.  7.5  15.6 

.  0 

.  0 

.  0 

77 . 0 

.  0 

1500.  50.  10000. 
0 


Data  input  file  for  RADCAL 

line  1:  number  of  homogeneous  elements,  n 

line  2:  pathlength  (m) ,  temperature  (K) ,  C02  (kPa) ,  H20,  CH4,  CO, 

02,  N2,  fv 

lines  3  through  n+1 :  same  as  line  2  for  the  rest  of  the  elements 
line  n+2 :  wall  temperature  (K) ,  minimum  wavenumber  (cm-1),  maximum 
wavenumber 

line  n+3:  0,  or  the  number  of  homogeneous  elements  in  the  next  case 
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APPENDIX  D.  TYPICAL  RADCAL  OUTPUT 


Radial  Profiles 


Partial  Pressures,  kPa 

J  dist, m  temp, K  C02  H20  CH4  CO 

1  .2850  1000.  7.500  15.600  .000  .000 

wall  400 . 


02  N2  FV 
.000  77.0  0.0 


Total  directional  radiated  energy  flux  =  .307139E+04  Watts/m-2/strad 


Spectral  Intensity  Distribution,  Watts/m-2/micron/strad 


micron 

intensity 

tau 

micron 

intensity 

tau 

1 . 005 

. 0000E+00 

1 . 0000 

6.061 

. 2463E+03 

.  8353 

1.010 

.  0000E+00 

1 . 0000 

6.154 

. 1477E+03 

.8975 

1 .015 

.  0000E+00 

1 . 0000 

6.250 

. 1126E+03 

.  9189 

1 . 020 

.  0000E+00 

1 . 0000 

6.349 

.  1516E+03 

.  8865 

1 . 026 

.  0000E+00 

1 . 0000 

6.452 

. 3674E+03 

.  7137 

1 . 031 

.  0000E+00 

1 . 0000 

6.557 

. 5144E+03 

.5824 

1 . 036 

. 0000E+00 

1 . 0000 

6.667 

.  4936E+03 

.5821 

1 . 042 

. 0000E+00 

1 . 0000 

6.780 

.  4175E+03 

.  6309 

1 . 047 

. 0000E+00 

1 . 0000 

6.897 

.  3891E+03 

.  6404 

1 . 053 

.  0000E+00 

1 . 0000 

7.018 

. 3447E+03 

.  6665 

1 . 058 

. 0000E+00 

1 . 0000 

7 . 143 

. 3104E+03 

.  6852 

1.064 

.  0000E+00 

1 . 0000 

7.273 

. 2909E+03 

.  6904 

1 . 070 

. 0000E+00 

1 . 0000 

7 .407 

. 2773E+03 

.  6897 

1 . 075 

.  0000E+00 

1 . 0000 

7 . 547 

. 2674E+03 

.  6850 

1.081 

.  2750E-03 

1 . 0000 

7 . 692 

.  2309E+03 

.  7133 

1 . 087 

.  5450E-03 

1 . 0000 

7 .843 

. 1981E+03 

.  7403 

1 .093 

. 5559E-02 

1 . 0000 

8 . 000 

. 1503E+03 

.7915 

1 .099 

. 7908E-02 

.  9999 

8 .163 

.  1068E+03 

.  8431 

1 . 105 

. 1017E-01 

.9999 

8 . 333 

.  7042E+02 

.8902 

1.111 

. 1929E-01 

.  9999 

8 . 511 

.  4206E+02 

.  9303 

1.117 

. 3087E-01 

.  9998 

8 . 696 

.  2218E+02 

.  9608 

1 . 124 

. 3884E-01 

.  9998 

8 . 889 

. 1643E+02 

.  9690 

1 . 130 

. 4722E-01 

.  9998 

9.091 

. 2560E+02 

.  9484 

1 . 136 

. 5675E-01 

.  9997 

9.132 

. 2  68 1E+02 

.  9452 

1 . 143 

. 6802E-01 

.  9997 

9 . 174 

. 2729E+02 

.  9434 

1.149 

. 7338E-01 

.  9997 

9.217 

. 2677E+02 

.  9437 

1 .  156 

. 9351E-01 

.  9996 

9.259 

. 2505E+02 

.9466 

1 .163 

. 1019E+00 

.  9996 

9.302 

. 2200E+02 

.  9524 

1 . 170 

.  1129E+00 

.  9995 

9.346 

. 17  68E+02 

.  9612 

1.176 

. 1128E+00 

.  9996 

9.390 

. 1214E+02 

.  9730 

1 .  183 

. 1404E+00 

.  9995 

9 . 434 

. 1413E+02 

.  9681 

1.190 

. 1451E+00 

.  9995 

9.479 

. 1813E+02 

.  9585 

1.198 

. 9713E-01 

.  9997 

9.524 

. 2072E+02 

.  9519 

1.205 

. 1105E+00 

.  9996 

9.569 

. 2199E+02 

.  9482 
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1 .212 

. 9183E-01 

.  9997 

9.615 

. 2207E+02 

.  9472 

1.220 

. 1101E+00 

.  9997 

9.662 

. 2119E+02 

.  9486 

1.227 

. 5030E-01 

.9999 

9.709 

. 1966E+02 

.  9515 

1.235 

. 4276E-01 

.9999 

9.756 

. 1780E+02 

.  9555 

1.242 

. 4418E-01 

.  9999 

9.804 

. 1598E+02 

.  9594 

1.250 

.  1951E-01 

1 . 0000 

9.852 

. 1440E+02 

.  9629 

1.258 

. 1816E-01 

1 . 0000 

9.901 

. 1324E+02 

.  9653 

1.266 

.  2371E-01 

.  9999 

9.950 

. 1255E+02 

.9666 

1 .274 

. 3800E-01 

.  9999 

10 . 000 

.  1225E+02 

.  9669 

1.282 

.  5689E-01 

.  9999 

10 . 050 

. 1224E+02 

.9664 

1.290 

. 7012E-01 

.  9999 

10 . 101 

. 1226E+02 

.  9658 

1.299 

.  3162E+00 

.  9994 

10 . 152 

. 1212E+02 

.  9657 

1 . 307 

. 3513E+00 

.  9993 

10 .204 

. 1167E+02 

.9664 

1.316 

. 8  4  66E+00 

.  9984 

10.256 

. 1080E+02 

.  9684 

1.325 

. 2019E+01 

.  9964 

10 . 309 

. 9539E+01 

.  9717 

1.333 

.  2680E+01 

.  9954 

10 . 363 

.  7912E+01 

.  9761 

1 . 342 

. 7485E+01 

.  9876 

10.417 

. 6633E+01 

.  9796 

1.351 

.  1995E+02 

.  9682 

10 . 471 

. 8154E+01 

.  9746 

1.361 

. 3008E+02 

.  9539 

10 . 526 

. 9307E+01 

.  9705 

1 . 370 

. 3041E+02 

.  9551 

10 . 582 

. 1006E+02 

.967  6 

1.379 

. 2758E+02 

.  9608 

10 . 638 

. 1035E+02 

.  9660 

1.389 

. 3325E+02 

.  9545 

10 . 695 

. 1024E+02 

.  9658 

1.399 

. 3583E+02 

.  9527 

10 . 753 

. 9810E+01 

.9667 

1 . 408 

. 3438E+02 

.  9563 

10 . 811 

. 9186E+01 

.  9683 

1.418 

. 3388E+02 

.  9585 

10 . 870 

. 8521E+01 

.  9701 

1.429 

. 30  91E+02 

.  9635 

10 . 929 

.  7865E+01 

.  9719 

1 . 439 

. 2601E+02 

.  9703 

10 . 989 

. 7284E+01 

.  9735 

1.449 

. 2227E+02 

.  9755 

11 . 050 

.  6814E+01 

.  9748 

1.460 

.  1790E+02 

.  9810 

11 .  Ill 

.  6460E+01 

.  9756 

1 . 471 

.  1061E+02 

.  9891 

11 . 173 

.  6270E+01 

.  9759 

1 . 481 

. 9673E+01 

.  9904 

11.236 

. 6159E+01 

.  9759 

1.493 

. 8050E+01 

.  9923 

11.299 

. 6104E+01 

.  9757 

1 . 504 

.  6925E+01 

.  9936 

11 . 364 

.  6029E+01 

.  9756 

1.515 

. 5957E+01 

.  9947 

11.429 

. 6057E+01 

.  9750 

1 . 527 

. 4139E+01 

.  9964 

11.494 

. 6112E+01 

.  9743 

1 . 538 

. 3552E+01 

.  9970 

11 . 561 

.  6162E+01 

.  9736 

1.550 

.  2128E+01 

.  9983 

11 . 628 

.  6207E+01 

.  9729 

1.563 

. 1522E+01 

.  9988 

11 . 696 

.  6248E+01 

.  9722 

1.575 

.  1290E+01 

.  9990 

11.765 

.  6283E+01 

.  9715 

1.587 

. 1014E+01 

.  9993 

11 . 834 

.  6376E+01 

.  9705 

1 . 600 

. 6685E+00 

.  9995 

11 . 905 

.  6463E+01 

.  9695 

1 . 613 

. 4852E+00 

.  9997 

11.976 

. 6542E+0 1 

.  9685 

1 . 626 

. 4411E+00 

.  9997 

12 .048 

. 7069E+01 

.  9653 

1 . 639 

.  4084E+00 

.  9997 

12 . 121 

. 7273E+01 

.  9636 

1 . 653 

.  38 66E+00 

.  9998 

12.195 

.  9147E+01 

.  9533 

1 . 667 

. 3688E+00 

.  9998 

12 .270 

.  8254E+01 

.  9570 

1 . 681 

.  4618E+00 

.  9997 

12.346 

. 8842E+01 

.  9530 

1 . 695 

. 5073E+00 

.  9997 

12 . 422 

.  9667E+01 

.  9475 

1.709 

. 6027E+00 

.  9997 

12.500 

.  1025E+02 

.  9432 

1 . 724 

.  1392E+01 

.  9993 

12.579 

. 1104E+02 

.  9376 

1 .739 

.  3258E+01 

.  9983 

12 . 658 

. 1197E+02 

.  9309 

1 .754 

. 9787E+01 

.  9950 

12.739 

. 1348E+02 

.  9205 

1 .770 

.  2302E+02 

.  9886 

12.821 

. 1497E+02 

.  9098 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


786 

. 4804E+02 

.  9769 

12 . 903 

. 1582E+02 

.  9026 

802 

. 8527E+02 

.  9601 

12 . 987 

. 17  67E+02 

.  8888 

818 

. 1165E+03 

.  9469 

13.072 

. 1978E+02 

.  8728 

835 

.  1231E+03 

.  9453 

13.158 

. 2243E+02 

.  8526 

852 

. 1330E+03 

.  9424 

13.245 

. 2550E+02 

.  8286 

869 

.  1840E+03 

.  9224 

13.333 

.  2466E+02 

.  8305 

887 

. 1082E+03 

.  9555 

13.423 

. 3142E+02 

.  7792 

905 

. 8995E+02 

.  9639 

13 . 514 

. 38  64E+02 

.  7222 

923 

. 1072E+03 

.  9580 

13.605 

. 3805E+02 

.  7201 

942 

. 1342E+03 

.  9487 

13.699 

. 3886E+02 

.  7074 

961 

. 1511E+03 

.  9435 

13.793 

.  47 69E+02 

.  6325 

980 

.  1217E+03 

.  9555 

13.889 

. 7088E+02 

.4409 

000 

. 1395E+03 

.  9501 

13.986 

. 6696E+02 

.  4591 

010 

. 91 92E+02 

.  9675 

14 . 085 

. 6327E+02 

.4766 

020 

. 8857E+02 

.  9690 

14 . 184 

. 6714E+02 

.4310 

030 

. 4371E+02 

.  9849 

14.286 

. 6781E+02 

.4113 

041 

.  4465E+02 

.  9847 

14 . 388 

. 7217E+02 

.3579 

051 

. 4845E+02 

.  9836 

14.493 

.  7992E+02 

.2711 

062 

. 3829E+02 

.  9872 

14 . 599 

. 8159E+02 

.2371 

073 

. 4135E+02 

.  9863 

14.706 

. 7775E+02 

.2545 

083 

. 1168E+02 

.  9962 

14.815 

.  8612E+02 

.  1531 

094 

.  4225E+01 

.  9986 

14 . 925 

. 9416E+02 

.  0500 

105 

. 1757E+01 

.  9994 

15.038 

. 90  91E+02 

.  0589 

116 

.  1581E+01 

.  9995 

15 . 152 

. 6921E+02 

.2646 

128 

. 1470E+01 

.  9995 

15.267 

. 6953E+02 

.2416 

139 

.  1308E+01 

.  9996 

15.385 

. 6904E+02 

.2267 

151 

.  1150E+01 

.  9996 

15.504 

. 5944E+02 

.  3162 

162 

.  1113E+01 

.  9997 

15.625 

. 5130E+02 

.3937 

174 

. 1056E+01 

.  9997 

15.748 

. 4844E+02 

.4116 

186 

.  1081E+01 

.  9997 

15.873 

. 4464E+02 

.4427 

198 

. 1131E+01 

.  9997 

16.000 

. 3964E+02 

.4913 

210 

.  1246E+01 

.  9996 

16.129 

. 3698E+02 

.5119 

222 

. 1342E+01 

.  9996 

16.260 

. 4188E+02 

.  4313 

235 

. 1390E+01 

.  9996 

16.393 

. 3814E+02 

.4672 

247 

.  2010E+01 

.  9994 

16.529 

. 3307E+02 

.5245 

260 

. 2084E+01 

.  9994 

16. 667 

. 3046E+02 

.5491 

273 

.  2685E+01 

.  9992 

16.807 

. 3175E+02 

.5159 

286 

. 3155E+01 

.  9991 

16.949 

. 2961E+02 

.5350 

299 

. 4134E+01 

.  9988 

17.094 

. 2789E+02 

.5486 

312 

. 4731E+01 

.  9987 

17.241 

.  2691E+02 

.5512 

326 

. 6020E+01 

.  9983 

17 .391 

.  2615E+02 

.5503 

339 

. 6801E+01 

.  9981 

17 . 544 

.  2582E+02 

.  5421 

353 

.  8293E+01 

.  9977 

17 . 699 

.  257 6E+02 

.5287 

367 

.  9228E+01 

.  9975 

17 .857 

.  2516E+02 

.5249 

381 

. 1570E+02 

.  9958 

18 . 018 

.  2490E+02 

.  5147 

395 

. 1983E+02 

.  9947 

18 . 182 

. 2457E+02 

.5054 

410 

. 3383E+02 

.  9910 

18.349 

. 2429E+02 

.4950 

424 

. 52  94E+02 

.  9860 

18 . 519 

.  2388E+02 

.4869 

439 

. 9285E+02 

.  9755 

18 . 692 

.  2350E+02 

.4780 

454 

. 1579E+03 

.  9586 

18 .868 

.  22  98E+02 

.4722 

469 

.  2717E+03 

.  9291 

19.048 

.  2256E+02 

.4641 

484 

. 4305E+03 

.  8883 

19.231 

. 2215E+02 

.  4556 

500 

. 6355E+03 

.  8359 

19.417 

. 2176E+02 

.4462 

73 


2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


516 

.  8206E+03 

.7891 

19.608 

. 2134E+02 

.  4374 

532 

.  9776E+03 

.7500 

19.802 

. 2090E+02 

.  4293 

548 

.  1067E+04 

.7282 

20 . 000 

. 2043E+02 

.  4217 

564 

. 9858E+03 

.7500 

20.202 

. 1998E+02 

.4136 

581 

. 1036E+04 

.7384 

20 . 408 

.  1952E+02 

.4060 

597 

.  7898E+03 

.  8012 

20 . 619 

.  1903E+02 

.3989 

614 

.  7146E+03 

.  8208 

20 . 833 

. 1853E+02 

.3923 

632 

. 5482E+03 

.8630 

21 . 053 

. 1803E+02 

.3861 

64  9 

. 4264E+03 

.8938 

21 .277 

. 1755E+02 

.3790 

667 

. 1657E+04 

.  5883 

21 . 505 

. 1706E+02 

.  3724 

685 

. 1849E+04 

.5419 

21 . 739 

.  1657E+02 

.3662 

703 

. 1975E+04 

.  5117 

21 . 978 

. 1607E+02 

.3603 

721 

. 1645E+04 

.5943 

22.222 

. 1557E+02 

.  3549 

740 

. 1340E+04 

.  6702 

22 .472 

.  1506E+02 

.  3499 

759 

.  1206E+04 

.7038 

22 . 727 

. 1456E+02 

.3452 

778 

.  1376E+04 

.  6625 

22 . 989 

. 1406E+02 

.3408 

797 

. 1440E+04 

.  6473 

23.256 

. 1356E+02 

.3366 

817 

. 1471E+04 

.  6400 

23.529 

. 1307E+02 

.  3326 

837 

. 1373E+04 

.6642 

23.810 

.  1258E+02 

.3291 

857 

. 1190E+04 

.7092 

24 .096 

. 1210E+02 

.3258 

878 

.  9948E+03 

.7570 

24 . 390 

. 1163E+02 

.3226 

899 

. 8855E+03 

.7837 

24 . 691 

. 1116E+02 

.3196 

920 

.  7522E+03 

.8163 

25.000 

. 1071E+02 

.3168 

941 

.64  4  6E+0  3 

.  8425 

25.316 

. 1025E+02 

.3150 

963 

.  5486E+03 

.8659 

25.641 

. 9802E+01 

.3134 

985 

.  4688E+03 

.  8853 

25 . 974 

. 9368E+01 

.3117 

008 

. 3842E+03 

.  9059 

26.316 

.  8945E+01 

.  3102 

030 

. 3267E+03 

.  9198 

2  6. 667 

. 8535E+01 

.3088 

053 

.  2719E+03 

.  9332 

27 . 027 

.  8122E+01 

.3086 

077 

.  2190E+03 

.  9461 

27 .397 

. 7724E+01 

.3085 

101 

. 1660E+03 

.  9590 

27 .778 

.  7339E+01 

.3084 

125 

.  1280E+03 

.  9683 

28 .169 

. 6967E+01 

.3083 

150 

. 1127E+03 

.  9720 

28 . 571 

. 6610E+01 

.  3082 

175 

. 9436E+02 

.  9765 

28 . 986 

. 6254E+01 

.3094 

200 

. 7517E+02 

.  9812 

29 . 412 

. 5913E+01 

.3105 

226 

. 6378E+02 

.  9840 

29.851 

. 5585E+01 

.3116 

252 

. 5180E+02 

.  9870 

30 . 303 

. 5271E+01 

.3127 

279 

. 4519E+02 

.  9886 

30 .769 

.  4971E+01 

.  3137 

306 

. 3939E+02 

.  9900 

31.250 

.  4673E+01 

.3162 

333 

. 3259E+02 

.  9917 

31.746 

.  4389E+01 

.  3186 

361 

. 2895E+02 

.  9926 

32.258 

. 4118E+01 

.3210 

390 

.  2597E+02 

.  9933 

32 . 787 

.  3860E+01 

.3234 

419 

. 2072E+02 

.  9946 

33.333 

. 3615E+01 

.3258 

448 

. 17  64E+02 

.  9954 

33.898 

. 3374E+01 

.3295 

478 

.  1464E+02 

.  9961 

34 . 483 

. 3145E+01 

.3333 

509 

.  1238E+02 

.  9967 

35.088 

. 2928E+01 

.  3371 

540 

.  1063E+02 

.  9972 

35.714 

.  2723E+01 

.3408 

571 

.  8686E+01 

.  9977 

36.364 

.  2529E+01 

.3445 

604 

. 6410E+01 

.  9983 

37 .037 

.  2338E+01 

.3505 

636 

. 4385E+01 

.  9988 

37.736 

. 2158E+01 

.3565 

670 

. 3877E+01 

.  9989 

38 .462 

.  1988E+01 

.3626 

704 

.  3277E+01 

.  9991 

39.216 

.  1829E+01 

.3687 

738 

.  2320E+01 

.  9993 

40 . 000 

.  1679E+01 

.  3749 
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3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 


774 

.  2023E+01 

.  9994 

40 .816 

. 1534E+01 

.3832 

810 

.  1604E+01 

.  9995 

41 . 667 

.  1398E+01 

.3918 

846 

. 1533E+01 

.  9996 

42.553 

. 1272E+01 

.4005 

883 

.  1524E+01 

.  9996 

43 . 478 

. 1154E+01 

.4095 

922 

. 1518E+01 

.  9995 

44 . 444 

. 1044E+01 

.4187 

960 

.  1460E+01 

.  9996 

45 . 455 

. 9412E+00 

.  4290 

000 

. 1480E+01 

.  9995 

46.512 

. 8458E+00 

.4395 

040 

.  1690E+01 

.  9995 

47.619 

.  7576E+00 

.  4505 

082 

. 1985E+01 

.  9994 

48 .780 

. 67  63E+00 

.4619 

124 

. 2304E+01 

.  9993 

50 . 000 

. 6015E+00 

.4737 

167 

.  2646E+01 

.  9991 

51.282 

. 5334E+00 

.4855 

211 

.  2953E+04 

.  0324 

52 . 632 

.  4711E+00 

.4978 

255 

.  2995E+04 

.  0031 

54 . 054 

. 4141E+00 

.5106 

301 

.  2954E+04 

.  0008 

55.556 

.  3623E+00 

.  5241 

348 

.  2905E+04 

.  0006 

57 . 143 

. 3152E+00 

.5383 

396 

.  2846E+04 

.  0039 

58 . 824 

. 2731E+00 

.5525 

444 

.  2625E+04 

.0648 

60.606 

.  2351E+00 

.5675 

494 

.  1643E+04 

.  4038 

62.500 

.  2009E+00 

.5835 

545 

.  57 66E+03 

.7868 

64.516 

. 1703E+00 

.  6006 

598 

. 1604E+03 

.  9395 

66. 667 

. 1430E+00 

.  6190 

651 

. 7032E+02 

.  9729 

68.966 

. 1200E+00 

.  6354 

706 

. 6562E+02 

.  9742 

71.429 

. 9956E-01 

.  6531 

762 

. 8193E+02 

.  9671 

74 . 074 

. 8162E-01 

.  6723 

819 

. 9943E+02 

.  9592 

76.923 

. 6593E-01 

.  6933 

878 

. 1257E+03 

.  9472 

80 . 000 

. 5227E-01 

.  7166 

938 

. 1503E+03 

.  9354 

83.333 

. 4202E-01 

.  7328 

000 

. 1674E+03 

.  9263 

86.957 

. 3324E-01 

.  7503 

063 

. 1852E+03 

.  9164 

90 . 909 

.  2578E-01 

.7695 

128 

. 1905E+03 

.  9118 

95.238 

.  1949E-01 

.7909 

195 

.  2603E+03 

.  8764 

100 . 000 

. 1424E-01 

.  8150 

263 

. 2817E+03 

.8625 

105.263 

. 1073E-01 

.  8295 

333 

.  3009E+03 

.  8491 

111.111 

. 7875E-02 

.  8452 

405 

. 3575E+03 

.8155 

117 . 647 

. 5586E-02 

.8625 

479 

. 3563E+03 

.8106 

125.000 

. 3781E-02 

.  8818 

556 

.  4000E+03 

.7808 

133.333 

. 2387E-02 

.  9038 

634 

. 4012E+03 

.7733 

142 . 857 

. 1644E-02 

.  9130 

714 

.  4043E+03 

.7642 

153.846 

. 1088E-02 

.  9228 

797 

.  4231E+03 

.7450 

166. 667 

. 6855E-03 

.  9333 

882 

. 3710E+03 

.7687 

181 .818 

. 4039E-03 

.  9445 

970 

.  2889E+03 

.8135 

200 . 000 

.  2162E-03 

.  9567 

The  effective  absorption  coef.  is  . 654587E+00/m 
The  Planck-mean  absorption  coef.  is  . 290665E+01/m 
The  wall-incident  mean  is  . 290665E+01/m 


Radial  Profiles 


Partial  Pressures,  kPa 
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J  dist, m  temp, K  C02  H20  CH4  CO  02  N2 

FV 

1  1.0000  1000.  7.500  15.600  .000  .000  .000  77.000 

. 0000E+00 

wall  0 . 

Total  directional  radiated  energy  flux  =  .125184E+04  Watts/m-2/strad 


Spectral  Intensity  Distribution,  Watts/m-2/micron/strad 


micron 

intensity 

tau 

micron 

intensity 

tau 

3.339 

. 1085E+03 

.  9723 

4 . 008 

. 5337E+01 

.  9984 

3.367 

. 9697E+02 

.  9751 

4.049 

. 6130E+01 

.  9981 

3.396 

. 8553E+02 

.  9779 

4 .090 

. 7178E+01 

.  9977 

3.425 

. 6933E+02 

.  9820 

4 . 132 

. 8308E+01 

.  9974 

3.454 

. 58  94E+02 

.  9846 

4 . 175 

. 2151E+04 

.  3038 

3 . 484 

. 4918E+02 

.  9870 

4.219 

. 3041E+04 

.  0006 

3.515 

. 4178E+02 

.  9889 

4.264 

. 2995E+04 

.  0000 

3.546 

. 3564E+02 

.  9905 

4 . 310 

. 2947E+04 

.  0000 

3.578 

. 2868E+02 

.  9923 

4 . 357 

. 2897E+04 

.  0000 

3.610 

. 2096E+02 

.  9943 

4 . 405 

. 2847E+04 

.  0000 

3.643 

. 1498E+02 

.  9959 

4 . 454 

. 27  93E+04 

.  0012 

3.676 

. 1314E+02 

.  9964 

4 . 505 

. 2490E+04 

.0932 

3.711 

. 1080E+02 

.  9970 

4 . 556 

. 1262E+04 

.5315 

3.745 

. 7920E+01 

.  9978 

4 . 608 

. 4181E+03 

.  8417 

3.781 

. 6793E+01 

.  9981 

4.662 

. 2144E+03 

.  9172 

3 . 817 

. 5574E+0 1 

.  9984 

4 . 717 

.  2115E+03 

.  9166 

3.854 

. 5368E+01 

.  9984 

4 . 773 

.  2523E+03 

.8983 

3.891 

. 5338E+01 

.  9984 

4 . 831 

.  2964E+03 

.  8778 

3.929 

. 5280E+01 

.  9984 

4 .890 

. 3505E+03 

.  8522 

3.968 

. 5131E+01 

.  9985 

4 . 950 

. 3933E+03 

.  8302 

The  effective  absorption  coef.  is  . 718888E-01/m 
The  Planck-mean  absorption  coef.  is  . 199213E+01/m 
The  wall-incident  mean  is  . 199213E+01/m 
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APPENDIX  E.  TYPICAL  TEP  OUTPUT 


INPUT  FILE :D: \PC_TEP\tep.dat 
OUPUT  FILE :D : \PC_TEP\tep. out 

THERMODYNAMIC  DATA  FILE : D : \PC_TEP\species\chno . dat 
TITLE  RPI/G02 
DATA 
$DATA 
ODE=l , 

NZONES=l , 

ASUB (1) =2, 

NASUB=1 , 

ASUP (1) =  6, 

NASUP=1 , 

$END 

REACTANTS 

298. 15F  .773 

298.150 

NAMELISTS 

$ODE 

RKT  =  .TRUE. 

PSIA=F, 

T (1) =1500, 

P  (1)=25, 

OFSKED (1) =3, 

OF=  T, 

$END 

OTITLE  RPI/G02 
ODATA 

$DATA 
ODE=l , 

NZONES=l , 

ASUB (1) =2, 

NASUB=1 , 

ASUP (1) =  6, 

NASUP=1 , 

$END 


C  1. 

0  2. 


H  1.9423 


100. 

100. 


-5430. L 
0.0G 


CALCULATE  ODE  AREA  RATIO  AND  PRESSURE  SCHEDULES  FOR  ZONE  1 


REACTANTS 


C  1.0000  H 

0.77300 

1.9423 

0.0000 

0.0000 

0.0000 

100.000000 

-5430.00 

L 

298.150 

O  2.0000 

0.0000 

0.0000 

0.0000 

0.0000 

100.000000 

0.00 

G 

298.150 

0.00000 

NAMELISTS 

$ODE 

RKT  =  .TRUE. 

PSIA=F, 

T (1) =1500, 

P ( 1 ) =25, 

OFSKED (1) =3, 

OF=  T, 

$END 

0SPECIES  BEING  CONSIDERED  IN  THIS  SYSTEM 


J  3/78 
9/65  C02 

C 

J12/67 

CH 

L  5/84 

CH4 

J  9/65 

CO 

J12/69 
6/88  JET- 

C2 

-A(G) 

J  3/61 

C2H2 

J12/69 

C3 

J12/69 

C4 

DR  9 

3/77  H2 

C12H26 

J  3/77 

H 

J12/70 

HCO 

J  9/78 

H02 

77 


J  3/79  H20 

J  3/77  02 

J  3/78  C (GR) 

PI  9  C12H2  6 (L) 

OOF  =  3.000000 

ENTHALPY 

(KG-MOL) (DEG  K) /KG 
OKG-ATOMS/KG 
C 
H 
0 

ENTHALPY  IN  BTU/LBM  : 
FROM  REACTANTS  : 
FROM  DELH (  )  : 

FROM  DELHI (  ) 
TOTAL  : 

1 

EXPANSION 


L  6/80  H202 

J  3/79  H20 (L) 


EFFECTIVE  FUEL 
HPP  (2) 

-0.19573987E+03 
BOP  (1,2) 

0 . 71587443E-01 
0 . 13904428E+00 
0.00000000E+00 


-174 . 9239 
0.0000 
0.0000 
-174 . 9239 

THEORETICAL  ROCKET 


J  3/77  O 

Est.  JP10(L) 


EFFECTIVE  OXIDANT 
HPP  (1) 

0.00000000E+00 
BOP  (1,1) 
0.00000000E+00 
0.00000000E+00 
0.62502339E-01 


J  6/77  OH 
L  6/88  JET-A(L) 


MIXTURE 

HSUBO 

-0 . 48934967E+02 
BO  ( I ) 

0 . 17896861E-01 
0 . 34761071E-01 
0 . 46876755E-01 


COMPOSITION  DURING 


ZONE  =  1 

PERFORMANCE  ASSUMING  EQUILIBRIUM 


OPC  =  367.4  PSIA 


FROM  AN  ASSIGNED  PRESSURE  AND  TEMPERATURE 

WT  FRACTION  ENTHALPY 


STATE  TEMP  DENSITY 

CHEMICAL  FORMULA  (SEE  NOTE)  CAL/MOL 

DEG  K  G/CC 

FUEL  C  1.00000  H  1.94230  1.00000  -5430.000 

L  298.15  0.7730 

OXIDANT  O  2.00000  1.00000  0.000 

G  298.15  0.0000 

0O/F=3 .0000E+00  PERCENT  FUEL=2 . 5000E+01  EQUIVALENCE  RATIO=l . 1343E+00  STOIC  MIXTURE  RATIO=3 . 4030E+00 
DENSITY=0 . 0000E+00 


1 

CHAMBER 

THROAT 

EXIT 

EXIT 

PC/P 

1.0000 

1.7704 

1.0598 

39.796 

P,  ATM 

25.00 

14.12 

23.59 

0.6282 

T,  DEG  K 

1500 

1365 

1486 

822 

H,  CAL/G 

-1845.1 

-1902 . 4 

-1851.2 

-2135.1 

S,  CAL/  (G)  (K) 

2.0645 

2.0645 

2.0645 

2.0645 

G,  CAL/GRAM 

-4941 . 9 

-4721 . 1 

-4918.6 

-3832.8 

U,  CAL/GRAM 

-1950.2 

-1998.1 

-1955.3 

-2192 . 6 

DEN  (G/LITER) 

5.76E+00 

3.57E+00 

5.48E+00 

2.65E-01 

M,  MOL  WT 

28.347 

28 . 347 

28.347 

28 . 427 

(DLV/DLP ) T 

-1.00000 

-1.00000 

-1.00000 

-1.00502 

(DLV/DLT ) P 

1.0000 

1.0000 

1.0000 

1.0833 

CP,  CAL/ (G) (K) 

0 . 4269 

0.4239 

0.4265 

0.5189 

CP  GAS(SF) 

0.4137 

0.4057 

0 .4129 

0.3650 

GAMMA  GAS (SF) 

1.2039 

1.2087 

1.2043 

1.2367 

GAMMA  (S) 

1.1963 

1.1980 

1.1965 

1 . 1806 

SON  VEL, M/SEC 

725.5 

692.6 

722.1 

532.9 

MU,  POISE 

5.24E-04 

4.89E-04 

5.20E-04 

3.32E-04 

K, ERG/S-CM-K 

1 . 41E+04 

1 .30E+04 

1 . 40E+04 

8.75E+03 

PRANDTL  NO 

0.64509 

0.63850 

0.64448 

0.57932 

MACH  NUMBER 

0.0000 

1.0000 

0.3124 

2.9246 

AE/AT 

1.0000 

2.0000 

6.0001 

CSTAR,  FT/SEC 

3358 

3358 

3358 

CF  VAC 

1.242 

1.673 

CF 

0.677 

1.522 

I VAC, LBF-S/LBM 

129.58 

174.65 

I,  LBF-SEC/LBM 

70.63 

158.91 

MOL  WT (MIX) 

28.347 

28 . 347 

28.347 

28 . 427 

MOLE  FRACTIONS 

CH4 

0.000000 

0.000000 

0.000000 

0.001415 

CO 

0 . 123704 

0.115658 

0 . 122957 

0.044893 

C02 

0.383614 

0.391660 

0.384361 

0.462446 

H 

0.000001 

0.000000 

0.000001 

0.000000 

H2 

0.054809 

0.062854 

0.055556 

0.128465 

H20 

0 . 437872 

0.429827 

0 . 437125 

0.362781 

MASS  FRACTIONS 
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CH4  0.000000  0.000000  0.000000  0.000799 

CO  0.122237  0.114286  0.121499  0.044235 

C02  0.595582  0.608074  0.596742  0.715947 

H  0.000000  0.000000  0.000000  0.000000 

H2  0.003898  0.004470  0.003951  0.009110 

H20  0.278283  0.273170  0.277808  0.229909 


OADDITIONAL  PRODUCTS 

ASSIGNED  CONDITIONS 

WHICH  WERE 

CONSIDERED  BUT 

WHOSE  MOLE 

FRACTIONS  WERE  LESS 

THAN  .0000005 

FOR  ALL 

C 

HCO 

CH 

H02 

C2 

C2H2 

C3 

C4 

JET-A(G) 

C12H26 

H202 

O 

OH 

02 

C  (GR) 

H20 (L) 

JP10 (L) 

JET-A(L) 

C12H26 (L) 


NOTE 

WEIGHT  FRACTION  OF  FUEL  IN  TOTAL  FUELS  AND  OF  OXIDANT  IN  TOTAL  OXIDANTS 
(SF)  STANDS  FOR  (SHIFTING  FROZEN) 

1 

0  FROZEN  TRANSPORT  PROPERTIES  CALCULATED  FROM  EQUILIBRIUM  CONCENTRATIONS 


CHAMBER 

THROAT 

EXIT 


MU 

(LBF-SEC/FT**2) 
1.09539496E-06 
1 . 021 99533E-06 
6 . 93548600E-07 


K 

(LBF/SEC-DEG  R) 
1 . 75734852E-02 
1 . 624  69596E-02 
1 . 09325694E-02 


VISCOSITY  EXPONENT  (OMEGA)  FOR  THE  FORM  MU=MUREF* (T/TREF)  1 


6 . 45086110E-01 
6 . 384  96220E-01 
5 . 7  9323411E-01 
^ OMEGA  IS  0.76373 


MUREF  FOR  INPUT  TO  BLM= 

SPECIES  CONSIDERED 
C 

C02 

H20 

02 


3 . 53317955E-05  LBM/ (FT-SEC) 

IN  TRANSPORT  PROPERTIES  CALCULATIONS 
CH  CH4  CO 

C2  C2H2  H2 

H202  O  OH 


ZONE  =  1 

THEORETICAL  ROCKET  PERFORMANCE  ASSUMING  FROZEN  COMPOSITION  DURING  EXPANSION 


0PC  =  367.4  PSIA 


FROM  AN  ASSIGNED  PRESSURE  AND  TEMPERATURE 

WT  FRACTION  ENTHALPY 


STATE  TEMP  DENSITY 

CHEMICAL  FORMULA  (SEE  NOTE)  CAL/MOL 

DEG  K  G/CC 

FUEL  C  1.00000  H  1.94230  1.00000  -5430.000 

L  298.15  0.7730 

OXIDANT  O  2.00000  1.00000  0.000 

G  298.15  0.0000 

0O/F=3 .OOOOE+OO  PERCENT  FUEL=2 . 5000E+01  EQUIVALENCE  RATIOl .  1343E+00  STOIC  MIXTURE  RATIO=3 . 4030E+00 
DENSITY=0 . 0000E+00 


1 

CHAMBER 

THROAT 

EXIT 

EXIT 

PC/P 

1.0000 

1.7776 

1.0602 

43.153 

P,  ATM 

25.00 

14.06 

23.58 

0.5793 

T,  DEG  K 

1500 

1359 

1485 

752 

H,  CAL/G 

-1845.1 

-1902.7 

-1851.2 

-2133.5 

S,  CAL/  (G)  (K) 

2.0645 

2.0645 

2.0645 

2.0645 

G,  CAL /GRAM 

-4941 . 9 

-4709.1 

-4917.5 

-3686.0 

U,  CAL /GRAM 

-1950.2 

-1997.9 

-1955.3 

-2186.2 

DEN  (G/LITER) 

5.76E+00 

3.57E+00 

5.48E+00 

2.66E-01 

M,  MOL  WT 

28.347 

28 . 347 

28.347 

28 . 347 

CP,  CAL/(G)(K) 

0.4137 

0.4047 

0 .4128 

0.3507 

GAMMA  (S) 

1.2039 

1.2093 

1.2044 

1.2496 

SON  VEL, M/SEC 

727.8 

694 . 4 

724.3 

525.0 

MACH  NUMBER 

0.0000 

1.0000 

0.3124 

2.9602 

AE/AT 

1.0000 

2.0000 

6.0002 

CSTAR,  FT/SEC 

3349 

3349 

3349 

CF  VAC 

1.243 

1.661 

CF 

0.680 

1.523 

I VAC, LBF-S/LBM 

129.36 

172.94 

I,  LBF-SEC/LBM 

70.81 

158.47 

MOLE  FRACTIONS 


CO  0.123704  C02  0.383614  H 

0.054809 


0.000001  H2 


79 


H20 


0 . 437872 


MASS  FRACTIONS 

CO  0.122237  C02  0.595582  H2  0.003898  H20 

0.278283 

0 ADDITIONAL  PRODUCTS  WHICH  WERE  CONSIDERED  BUT  WHOSE  MOLE  FRACTIONS  WERE  LESS  THAN  .0000005  FOR  ALL 
ASSIGNED  CONDITIONS 


c 

HCO 

CH 

H02 

C2 

C2H2 

C3 

C4 

JET-A(G) 

C12H26 

H202 

C12H26 (L) 

O 

OH 

02 

C  (GR) 

H20 (L) 

JP10 (L) 

JET-A(L) 

NOTE 

WEIGHT  FRACTION  OF  FUEL  IN  TOTAL  FUELS  AND  OF  OXIDANT  IN  TOTAL  OXIDANTS 

1 


****************************************************************************************************** 

****■*■*■*■*■*■*■*■*■*■*■*■*■*■* 

CALCULATE  ODE  AREA  RATIO  AND  PRESSURE  SCHEDULES  FOR  ZONE  1 

****************************************************************************************************** 

******■*■*■*■*■*■*■*■*■*■*■*■* 

***  EOF  ENCOUNTERED  IN  READING  ODE  REACTANTS  DATA  *** 

*********CPU (sEC)  =  0.0  ********** 
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